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We introduce a method that allows the evaluation of general expressions for the spectral functions 
of the one- dimensional Hubbard model for all values of the on-site electronic repulsion U. The 
spectral weights are expressed in terms of pseudofermion operators such that the spectral functions 
can be written as a convolution of pseudofermion dynamical correlation functions. Our results are 
valid for all finite energy and momentum values and are used elsewhere in the study of the unusual 
finite-energy properties of quasi-one-dimensional compounds and the new quantum systems of ultra- 
cold fermionic atoms on an optical lattice. 
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systems; heavy fermions 



I. INTRODUCTION 



The main goal of this paper is to provide a general method for the evaluation of matrix elements of one-, two- 
electron, or A/"-electron operators such that J\f is finite, between the ground state and excited energy eigenstates of 
the one-dimensional (ID) Hubbard model. Our results correspond to an important part of the derivation of the 
one-electron and two-electron spectral- weight distributions used in Refs. PJQ m the study of the phase diagram and 
unusual one-electron spectral properties of quasi-lD compounds. Indeed, the matrix-element and general spectral- 
function expressions derived here are used in Ref . Q in the evaluation of closed-form expressions for the one-electron 
and two-electron spectral-weight distributions of the model metallic phase. The ID Hubbard model is one of the few 
realistic models for which one can exactly calculate all the energy eigenstates and their energies 0, ■ In addition 
to the applications to the study of the unusual properties of the quasi-lD compounds presented in Refs. 0,0, our 
results are also of interest for the understanding of the spectral properties of the new quantum systems described by 
ultra-cold fermionic atoms in optical lattices with on-site repulsion |(| . 

The electron - rotated-electron unitary transformation used in Ref. @ for all values of the on-site repulsion 
U and the pseudofermion scattering theory introduced in Ref. play a central role in the construction of the 
pseudofermion dynamical theory introduced here. For finite values of U very little is known about the finite-energy 
spectral properties of the model. This is in contrast to simpler models [lOj. Unfortunately, combination of the 
model Bethe-ansatz solution 0, [j| with bosonization, conformal-field theory, or g-ology and Renormalization Group 
[lliri2| only allows the derivation of low-energy correlation-function expressions. In the limit of infinite U the spectral 
functions can be evaluated by the method presented in Ref. 13] and there are recent numerical results for finite values 
of U | 14| , but it is difficult to extract from them information about the microscopic processes that control the unusual 
spectral properties of the model. 

The paper is organized as follows: In Sec. [H] we introduce the model and the spectral-function problem and 
summarize the pseudofermion operational description used in our study. In Sec. IIIII we write the general spectral 
functions in terms of rotated-electron operators. The description of the rotated-electron elementary processes in terms 
of pseudofermion operators is the problem addressed in Sec. I1VI In Sec. ]V\ we express the spectral functions as a 
convolution of pseudofermion spectral functions and study the pseudofermion determinants involved in the expressions 
of these functions. Finally, the concluding remarks are presented in Sec. I VI I 

II. THE MODEL, THE SPECTRAL FUNCTIONS, AND THE PSEUDOFERMION DESCRIPTION 

In a chemical potential fi and magnetic field H the ID Hubbard Hamiltonian can be written as, 

H = H so[i) + V* S a z ■ Hso(A) =-t£E E 4 c c j+s, ° + U E T - V2] [ h jA - 1/2] . (1) 

a=c,s j = l<7=|,i5=-l, + l j=l 
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Here the operators c] a and Cj t „ are the spin-projection a electronic creation and annihilation operators at site j and 
fij. a = cjj a Cj. cr where j = 1,2, ...,JV . The number of lattice sites iV a is even and very large. We consider periodic 
boundary conditions. In the first expression of Eq. QJ, /x c = 2/x, /U s = 2/xoiT, /j,q is the Bohr magneton, and the 
diagonal generators of the fy-spin and spin SU(2) algebras S£ and S|, respectively, are given in Eq. (2) of Ref. 

The Hamiltonian Hso{a) of Eq. (JTJ commutes with the six generators of these two algebras, their off-diagonal 
generators being given in Eqs. (7) and (8), respectively, of Ref. The electron number operator reads N — N a , 
where N a = Ylf—x hj.a- F° r simplicity, we use units such that the Planck constant and electronic lattice constant are 
one. The model Q describes iVj spin-up electrons and spin-down electrons in a chain of N a sites, whose length 
in the units used here reads L — N a . We introduce the Fermi momenta which, in the thermodynamic limit L — > oo, 
are given by ±&f<t — ±7T"<t and ±kp = ±[kp^ + fcpjJ/2 = ±7m/2, where n a = N a /L and n = N/L. The electronic 
density can be written as n = n.| + n j and the spin density is given by m = n.| — n.| . We denote the ?7-spin value 
r] and projection rj z (and spin value S and projection S z ) of an energy eigenstate by S c and S 1 ,? (and S s and 5*1), 
respectively. The momentum operator reads, 



+7T 

dkN a {k)k, (2) 



<r=T,l fc <7=T,f w 

and commutes with the Hamiltonians given in Eq. (I}. The spin-projection tr momentum distribution operator 
appearing in Eq. (J2J is given by N a (k) = c\ ff Cfc jCr . Here CT and Cfe i<T are the spin-projection <r electron creation 
and annihilation operators of momentum k. These operators are related to the above local operators as follows, 

1 N a N a 

t = _L_ \^ A ■ ,,. . - _L_ g e"^ a c,- a . (3) 



The Bethe-ansatz solvability of the ID Hubbard model is restricted to the Hilbert subspace spanned by regular 
states, i.e. the lowest- weight states (LWSs) of the ?7-spin and spin algebras such that S a — —S", where a = c, s |. 
For simplicity, in this paper we restrict our considerations to values of the electronic density n and spin density m 
such that < n < 1 and < m < n, respectively. Often our expressions are different for the n = 1 Mott-Hubbard 
insulator phase and < n < 1 metallic phase (and for m — zero spin density and < m < n finite spin densities). 

The main aim of this paper is the evaluation of expressions for finite-w A^-electron spectral functions Bj^(k, ui), 
such that I = ±1, of the general form, 



(k,u) = ^ \{f\d l M (k)\GS)\ 2 s(u>-l[E f -E GS ]); lu;>0; I = ±1 , (4) 



where the operators in the matrix elements are such that, 

6+1 {k) = 6Uk) ; 6 M \k) = 6 N (k) . (5) 

Here the / summation runs over the excited energy eigenstates, the energies Ef correspond to these states, Eqs is 
the ground-state energy, and we use a momentum extended scheme such that k £ (-co, +oo). The operators OL-(fe) 
and Oj^(k) carry momentum k and are denoted in Eq. (JSJ by Ov(fc) where I = +1 and I = — 1, respectively. They 
are related to the local operators 6\f ■ = Ofjj and 6a", j = OJ^ 1 ■, respectively, by a Fourier transform. 
The local operators O 1 ^ j can be written as a product of 

l c , I,=±l 

local electronic creation and annihilation operators. Here J\fj l is the number of local electronic creation and 
annihilation operators of the operator Ojy • for l c = — 1 and l c = +1, respectively, and with spin down and 
spin up for l s = —I and l s = +1, respectively. It is assumed that the ratio Af/N a vanishes in the thermody- 
namic limit. Note that, by construction, Af + l and Af[~ t are such that Af+ l = NZi j j N-x i > A^'^, and 
NZ\ i < ■A/'+i ; • F° r A/" > 1 the operator Ojy ^ has a well defined local structure involving the Af-x,i 3 elec- 
tronic creation operators of spin projection l s /2, and M l +1 ; electronic annihilation operators of spin projection 
l s /2 located in neighboring lattice sites. The more usual cases for the description of experimental studies corre- 
spond to N = 1 and TV = 2. Examples of TV-electron operators 0_\f(k) = Oj/(k) are the one-electron operator 
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0\{k) = Ck,a (measured in the angle-resolved photoelectron spectroscopy), the spin-projection a density operator 
0% sd (k) = -j==J2k> c l+k>,a c k',cr> the transverse spin-density operator d| du, (/c) = ^= J2k> c l+k>. i c k'A, the on-site 
s-wave singlet superconductivity operator 02 SS (k) — -^=J2k ,c k-k',TCk\i, and the spin-projection <r triplet super- 
conductivity operator 02* s (fc) = -j^=J2k' cos(k')ck-k\crCk',a- The corresponding local operators Ojjj = Omj are 
6ij = Cj l(T , O^f = c^c-j.v, O^ff = ct^c,-,;, 0|*j = Cj, T Cj,|, and 00 = Cj )0 -Cj+i l0 -, respectively. The charge 
density operator (measured in density-density electron energy loss spectroscopy and inelastic X-ray scattering) is 
written in terms of the above spin-up and spin-down density operators. The operators O^(fc) of physical interest, 
correspond in general to operators O 1 ^ j whose Af elementary electronic operators create or annihilate electrons in a 
compact domain of lattice sites. For instance, if O^j = c ] j. c l+i i and thus O^* = Cj+j. f j, the interesting cases 
correspond to i = (on-site s-wave singlet superconductivity) and i = 1 (extended s-wave singlet superconductivity). 

The k dependence of the spectral functions Q can be transferred from the A/"-electron operators O l A f(k) to the 
excited energy eigenstates as follows, 

B^(k 1 Lo)^Y, N -\(f\ d ^,o\ GS )\ 2s (^- l i E f- E Gs])Sk,iik^kas]^ iw>0; i = ±l. (7) 

Here, Ojy is the J = local operator 0\f ■ considered above, fc/ is the momentum of the excited energy eigenstates, 
and kcs denotes the ground-state momentum. In this expression, we have chosen j = for the local operator Oj^ • . 
Due to translational invariance, the value of the functions 10 is independent of this special choice. 

Let us summarize the basic information about the holon, spinon, pse udop article, and pseudofermion descriptions 
needed for our studies. (For further information, see Refs. 0, H, 0, Il8|.) These studies involve the electron - 
rotated-electron unitary transformation, such that rotated-electron double occupancy is a good quantum number for 
all U/t values Q. As the Fermi- liquid quasiparticles, the rotated electrons have the same charge and spin as the 
electrons, but refer to all energies and reorganize in terms of [N a — N c ] 77-spin 1/2 holons, N c spin 1/2 spinons, and 
N c spinless and ?y-spinless c pseudoparticles, where N c is the number of rotated-electron singly occupied sites 0]. 
We use the notation ±1/2 holons and ±1/2 spinons, which refers to the 7/-spin and spin projections, respectively. 
The ±1/2 holons of charge ±2e correspond to rotated-electron unoccupied (+) and doubly-occupied (— ) sites. The 
complex behavior occurs for the spin-projection er-rotated electrons occupying singly occupied sites: their spin degrees 
of freedom originate chargeless spin-projection a spinons, whereas their charge part gives rise to 77-spinless and spinless 
c pseudoparticles of charge — e. 

Based on symmetry considerations, we can classify the ±1/2 holons and ±1/2 spinons into two classes: those which 
remain invariant under the electron - rotated-electron unitary transformation, and those which do not. The former 
are called independent ±1/2 holons and independent ±1/2 spinons. For instance, the ±1/2 Yang holons and ±1/2 
HL spinons 0,H,0)0| with numbers reading L c ±1 / 2 — [S c =F S"] and £ Sl ±i/2 = [S a T S%], respectively, belong to 
the former group of holons and spinons. The latter are part of ?7-spin-zero 2^-holon composite cv pseudoparticles and 
spin-zero 2j/-spinon composite sv pseudoparticles, respectively, where v = 1,2,... is the number of +1/2 and —1/2 
holon or +1/2 and —1/2 spinon pairs. Thus, the total number of ±1/2 holons (a = c) and ±1/2 spinons (a = s) 
reads M a ±1 / 2 = L a ,±i/2 + S^Li v ^avi where N av denotes the number of composite av pseudoparticles. The total 
number of holons (a — c) and spinons (a = s) is then given by M a = L a + 2 X)^°=i v N a v where L a = 2S a denotes 
the total number of Yang holons (a = c) and HL spinons (a — s). These numbers are such that M c =_\N a — N c ] and 
M s = N c . The pseudoparticles can be defined in terms of bare-momentum or spatial coordinates |17| . In addition 
to the Yang holons and HL spinons, also the holons and spinons associated with av ^ cO, si pseudoparticles of 
limiting bare-momentum values ±q av are independent holons and spinons. {q av is given in Eq. (B.14) of Ref. |8|].) 
Indeed, the invariance under the electron - rotated-electron unitary transformation of such cv pseudoparticles (and sv 
pseudoparticles) implies that they separate into 2v independent holons (and 2v independent spinons) and a cv (and 
sv) FP scattering center (These centers are defined in Ref. 0.) The emergence of the exotic quantum phases 
of matter considered in our study involves a second unitary transformation, which maps the c pseudoparticles (and 
composite av pseudoparticles) onto c pseudofermions (and composite av pseudofermions) . Such a transformation 
introduces shifts of order 1/L in the pseudoparticle discrete momentum values and leaves all other pseudoparticle 
properties invariant. Here we use the designation cO pseudoparticle and pseudofermion for the c pseudoparticle and 
pseudofermion, respectively. Thus, the cv and sv branches are such that v = 0,1, 2, ... and v — 1, 2, respectively. 

The local av pseudofermion creation (and annihilation) operator /J av (and f Xj , av ) creates (and annihilates) 

a av pseudofermion at the effective av lattice site of spatial coordinate xj = j a^. Here j = 1, 2, N* v and 
= L/N*^ = Na/N*^ is the effective av lattice constant given in Eq. (55) of Ref. 01 i n units of the electronic 
lattice constant. The general expression of the number of effective av lattice sites iV*^ is given in Eq. (B6) of 
Ref. 8], where the number of av pseudofermion holes N^ v is provided in Eq. (B.ll) of the same reference. (The 
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number of pseudofermion and pscudofermion holes equals that of the corresponding pseudoparticle and pseudoparticle 
holes The operator j\. av (and f Xjtav ) is related to the operator au (and fq j>a u)> which refers to av 

pseudofermions of canonical- momentum q~j , by a Fourier transform. The discrete canonical- momentum values of the 
av pseudofermions have a functional character and read , 

Qi=1j + Qt(<lj)/L=[27r/L}I^ + Qt(q j )/L; j = 1, 2, .., N* av , (8) 

where qj = [2irI° LV \/L is the bare-momentum carried by the av pseudoparticles. Here I" v are the actual quantum 
numbers provided by the Bethe-ansatz solution [8j. Although the av pseudoparticles carry bare- momentum qj, one 
can also label the corresponding av pseudofermions by such a bare-momentum. When one refers to the pseudofermion 
bare- momentum qj, one means that qj is the bare-momentum value that corresponds to the canonical momentum 
= qj + Q% v (qj)/L. Here and in Eq. (JBJ) Qtt V {qj)/2 i s a av pseudofermion overall scattering phase shift given by 

Qtfe)/2 = 7T^ J2 V Q vfe,9 3 0AJV aV fe); j = l,2,...,N* v , (9) 

a'v' j' = l 

where AN au (qj) = AAf a ^((jj) is the distribution function deviation AN al/ (qj) — N alJ (qj) — N®„(qj). The canonical- 
momentum distribution function J\f a v{q~j) (and bare- momentum distribution function N ai/ (qj)) is given by M a v (qj ) = 1 
and N au (qj) = (and N av (qj) — 1 and N n[/ (qj) = 0) for pseudofermions and pseudofermion holes (and pseudoparticles 
and pseudoparticle holes), respectively [3j. The ground-state densely- packed bare-momentum distribution function 
Nav(lj) i s defined in Eqs. (C.1)-(C.3) of Ref. 0. The av ^ cO, si pseudofermion limiting canonical-momentum 
values play an important role in the theory and read, 

q° v = lav + Qtuilau) / L ; ai/^cO, si, (10) 

where q av is the ground-state limiting bare-momentum value given in Eqs. (C.13) and (C.14) of Ref. |3] and q av the 
excited-energy-eigenstate limiting bare-momentum value provided in Eq. (B.14) of the same reference. In contrast 
to the av pseudoparticles, the av pseudofermions have no residual-interaction energy terms 0. Instead, under 
the ground-state - excited-energy-eigenstate transitions the av pseudofermions and av pseudofermion holes undergo 
elementary scattering events with the a'v' pseudofermions and a'v' pseudofermion holes created in these transitions 
0. This leads to the elementary two-pseudofermion phase shifts 7r $ Q „. a 'v'(qj ', q'j) on the right-hand side of the 
overall scattering phase shift J^J), which are defined by a set of integral equations |3J,|9|. The overall av pseudofermion 
or hole phase shift, 

Qa,fe)/2 = QL/2 + QLfe)/2, (11) 

plays an important role in the pseudofermion theory [10. Here Q°„/2 can have the values Q° av /2 = 0, ±tt/2 0,0. 

The pseudofermion description refers to a Hilbert subspacc called pseudofermion subspace (PS) where the Af- 
electron excitations are contained 0. The PS is spanned by the initial ground state and the excited energy 
eigenstates generated from it by the following types of processes (A)-(C), which are defined in more detail in Ref. 0]: 
(A) finite-energy and finite-momentum elementary cO and si pseudofermion processes plus creation of av ^ cO, si 
pseudofermions with bare-momentum values q ^ ±g QJ y, (B) zero-energy and finite- momentum processes that change 
the number of cO and si pseudofermions at their Fermi points, which for the ground state and L — > oo read, 

q° Fc0 = 2k F ; q° Fsl = k F[ , (12) 

plus creation of independent —1/2 holons and/or —1/2 spinons, and (C) low-energy and small-momentum elementary 
cO and si pseudofermion particle-hole processes in the vicinity of the Fermi points. The PS contains subspaces of 
several CPHS ensemble subspaces. (Here CPHS stands for cO pseudofermion, holon, and spinon.) A CPHS ensemble 
subspace is spanned by all energy eigenstates with fixed values for the —1/2 Yang holon number L c . _i/2, —1/2 HL 
spinon number £3,-1/2, cO pseudofermion number iV c o, and for the sets of composite av pseudofermion numbers 
{N cy } and {A^^} corresponding to the v = 1,2, ... branches. 

The pseudofermion bare-momentum dependent energy dispersions e c o(q), e s i(q), e cu (q) = 2vji + e^ v (q) for v > 0, 
and t S v(q) = 2v[1qH + e® u (q) for v > 1, where \i = /x(n) and H = H(m) correspond to the density and magnetization 
curves, respectively, are denned and studied in Refs. 000. Such energy dispersions play a crucial role in the 
expressions of the A/"-electron spectral functions. For m = 0, the energy 2\i is an increasing function of U and a 
decreasing function of the density n with the following limiting values, 

2/i = 4icos(7ma/2) , U/t -> ; U + it cos(ir na) , U/t — >• 00 ; U + At , n^0; E MH , n ->• 1 , (13) 
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where Emh is the half- filling Mott-Hubbard gap [4(. 

The evaluation of matrix elements between energy eigenstates considered in Sec. V involves pseudofermion operators 
fq. av an d fq\ av such that the canonical momentum values q and q' = q' correspond to an excited-energy-eigenstate and 
initial ground-state CPHS ensemble subspaces, respectively. In that case the unusual pseudofermion anticommutation 
relations read |2l|9j, 

i sinf Q alJ (q)/2 ) 

/ft f, , ,\-S c -ifg-g'V2 iQ„„(c,V2 V 6 " w/ J (u) 

N* v sin([g - q']/2) 

an d {fq t avi fq' t a',!/'} = {/g, af; /g', a', v'} = 0. 



III. SPECTRAL FUNCTIONS IN TERMS OF ROTATED-ELECTRON OPERATORS 

Here we express the general A/"-electron spectral functions Q in terms of rotated-electron creation and annihilation 
operators and evaluate the spectral- weight contributions from the Yang holons and HL spinons. Our first goal is 
identifying the set of CPHS ensemble subspaces which are spanned by the excited energy eigenstates generated by 
application onto the initial ground state of the operator O 1 ^ of Eq. Q . For clarification of this problem, we must 
find the set of deviation numbers AY c0 , AY sl , {AL a _ 1 / 2 }, and {AN a „} for av ^ cO.sl that are generated by 
application onto the ground state of that operator. According to the results of Refs. jjj for the ground state 
M c _i/ 2 = L a _i/ 2 = N al/ = for the av ^ cO, si branches and thus AM C .-1/2 = M Cj _i/2, AL Q , _ 1 y 2 = i a!) -i/2) 
and AN aL , = N av for the latter branches. 

First, we note that the values of the +1/2 holon and +1/2 spinon number deviations are such that, 

AM c , +1/2 = -AN cQ - M C) _i /2 ; AM S:+1/2 = AN cQ ~ AM s ,_ 1/2 , (15) 

and thus are dependent on the values of the —1/2 holon and —1/2 spinon numbers and cO pseudofermion number 
deviations. Also the occupancy configurations of the —1/2 holons and —1/2 spinons determine those of the +1/2 holons 
and +1/2 spinons. Indeed, the —1/2 holons and +1/2 holons correspond to the rotated-electron doubly-occupied sites 
and unoccupied sites, respectively, of a charge sequence 0. The point is that the spatial position of the unoccupied 
sites corresponds to the sites left over by the rotated-electron doubly occupied sites of a charge sequence. The same 
applies to the —1/2 spinons and +1/2 spinons, provided that we replace the rotated-electron doubly-occupied sites 
and unoccupied sites by sites singly occupied by spin-down and spin-up rotated electrons, respectively, and the charge 
sequence by the spin sequence [T3| - Moreover, the values of the corresponding +1/2 Yang holon and +1/2 HL spinon 
number deviations read, 

00 

AL C , +1/2 = -AN cQ -2J2"N C<1/ - L Ct _ 1/a , (16) 
i/=i 

and 

00 

AL S , +1/2 = AiV c0 - 2AiV sl - 2 ^ v N s , „ - L s , _ 1/a , (17) 

v=2 

respectively, and thus are not independent. One does not need to provide these values in order to specify a CPHS 
ensemble subspace. Therefore, often we do not consider in the expressions below the values of the holon numbers 
Mc.+i/2 an d I c .+i/2 and of the spinon numbers M s +1 / 2 and L s +1 / 2 - 

The values of the deviations AJVj and AN^ specific to a given A/"-electron operator, lead to sum rules for the values 
of the number deviations of pscudofcrmions, —1/2 Yang holons, and —1/2 HL spinons as follows, 

00 

AY = AY c0 + 2L c _ 1 /2 + 2^ v N cv , (18) 

and 



00 

A(Y X - Y T ) = 2AY sl - AY c0 + 2L s> _ 1/2 + 2^"^ 

v=2 



(19) 
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Only transitions to excited energy eigenstates associated with deviations obeying the sum rules l|18[) and l|19|l are 
permitted. The same deviations are associated with sum rules obeyed by the numbers A/"/ ; of Eq. JBJ) for the 

operator Oj^f • appearing in the general spectral- function expressions of Eq. J2J ■ Such sum rules read, 

AN = E HX,i.i A(N l -N ] )= J2 (20) 

l e ,l s =±l z c ,; s =±i 

Furthermore, it is straightforward to show that the following selection rule is valid for initial ground states corre- 
sponding to the density values considered in this paper: The values of the numbers of —1/2 Yang holons and —1/2 
HL spinons generated by application onto the ground state of the TV-electron operator O 1 ^ ■ , Eq. Q , are restricted 
to the following ranges, 

L C) _ 1/2 = 0,1,2,..., M-his> ^,-1/2 = 0,1,2,..., J2 5 I = ±1 » ( 21 ) 

ls=±l lc,ls=±l 

respectively. Here the numbers TV/ ; are those of Eq. @J specific to that operator. 

Further selection rules in terms of the rotated-electron expressions for the operator O 1 ^ of the general spectral- 
function |0 are given in the ensuing section. 

Let us label the excited energy eigenstates of the state summations of the general A/"-electron spectral function Q 
according to their CPHS ensemble subspace. (We recall that all excited energy eigenstates belonging to a given CPHS 
ensemble subspace have the same values for the set of deviation numbers A7V c0 and A7Y s i and numbers {L a _ 1 / 2 } for 
a = c, s, and {N av } for the av ^ cO, si branches.) This procedure leads to the following expression for the spectral 
function J7J), 

B l „(k,w)= J2 N ^f^ C \ & ^,o\GS)\ 2 6(^-lAE cph ^S k j Akcphs ; leu > , (22) 

{AN W }, -i/a} / 

where I = ±1 and the summation X){atv } {l t 2 } runs over the CPHS ensemble subspaces whose deviation values 
obey the sum rules (|18f) and (|19fl and selection rules l|21|l . Moreover, the summation J2f runs over the excited energy 
eigenstates |/; C) of a given CPHS ensemble subspace, AE cp h s is the excitation energy, and Ak cp h s the corresponding 
excitation momentum. A general energy eigenstate |/) with finite values for the numbers L c> _i/a and/or L s _i/ 2l 
can be expressed as follows, 

|/)= J] (S « )L l: 1/2 \f.L). (23) 



L, 



, -1/2 



Here, 

C a = <5l q ,_ 1/2 ,o+ ]J l[L a + l-l]\ L a ^ 1/2 <L a = 2S a , (24) 

i=i 

and the 77-spin flip Yang holon (a — c) and spin flip HL spinon (a = s) operators are the off-diagonal generators 
of the corresponding 5/7(2) algebras given in Eqs. (7) and (8), respectively, of Ref. @- These operators remain 
invariant under the electron - rotated-electron unitary transformation and thus have the same expression in terms 
of electronic and rotated-electron creation and annihilation operators. Moreover, in Eq. (|23|) \f.L) is the LWS that 
corresponds to the state |/). For a state \ f;C) belonging to a given CPHS the corresponding LWS is denoted by 
\f.L ; C). However, note that a non-LWS |/ ; C) and the corresponding LWS \f.L ; C) belong to different CPHSs, once 
they correspond to different values of the numbers £ c , -1/2 and/or L s _ l / 2 . 

It is useful to reexpress the spectral-function expression i|22[) in terms of matrix elements between regular states 
only. The ground state is a LWS of both the 77-spin and spin SU(2) algebras and thus has the following property, 

S a \GS)=0; a = c,s. (25) 

Let us introduce the operators 0jy • and 8^ k such that, 

(f.L;C\ J] * (S Q )^.-v 2 o^J G s) = - *-\(f.L ; C\&„JGS), 

ol—c, s v a— c, s v 

(f.L;C\ [] -^={S a ) L ^-^& M {k)\GS) = [[] -^=\(f-L;C\e'^ k \GS); l = ±l. (26) 
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By suitable use of Eq. I|25|) , it is straightforward to show that the operators 8jy ■ and 8^- k are given by the following 
commutators, 

©k, = [ n ^«) £ — i/a . 



e 



JV.fc 



[ II (S a ) La -- 1/2 , OV.fc] . L c , _i/2 and/or L 5 ,_ 1/2 >0; f = ±1 , (27) 



or by, 



Thus, 



®jv,j = o 1 at,3 ®jv, fc — Oa/ - , fc ) ^ c ,-i/2 = £ s ,-i/2 = 0; i = ±l. (28) 

B l N {k,u>) = J2 ( II r) J2 N ^f- L ' C \^o\GS)\ 2 

{AJf„,},{£ a ,_ 1/a } «=c,s Q / 

x <s(c- A£ cp/ls ) <5fe,Afc cph3 ; l = ±l. (29) 



Note that when the operator 0jy is given by Eq. (|28[) one has that \f.L; C) = \ f: C) in Eq. i|29|) . If the commutator 
iria=c s (S a ) Lc, '~ 1/2 , O 1 ^- ■] of Eq. 1(27)1 vanishes, then the excitation generated by application of the corresponding 
operator Ojy ■ onto the ground state has no overlap with the excited energy eigenstate <|23[) . 

Similarly to O 1 ^ ■, the corresponding operator 0^- . can be written as a A/"-electron operator. Let the numbers 
A/J j of Eq. © refer to the operator Ojv j of the spectral-function expression of Eq. (J7J. Then we call N{ i s the 
corresponding electronic numbers of the operator 0^ ■ of the spectral- function expression of Eq. (|29|l . The values 
are such that = £j , =±i-#? , = J2i c ,i t =±i-^Li s and > 



E A7 ±Ua= E ^±i,l. ±2£ o.-Va; E *.,±'.K,i.= E *.,±i.<,i. : F2i,,_i /a ; i = ±l. (30) 

These relations provide information about the numbers of electronic creation and annihilation operators of the Af- 
electron operator • expression relative to the corresponding numbers of the O 1 ^ j expression. While the number 
of electronic creation (and annihilation) operators decreases (and increases) by 2L c _i/2, the number of electronic 
spin-down creation and spin-up annihilation (and spin-down annihilation and spin-up creation) operators decreases 
(and increases) by 2L s _i/ 2 - 

We note that the numbers Af[ a 1b of Eq. © for the operator 0\f 3 on the right-hand-side of Eq. \27\ obey the sum 

rules H2Ufl . Thus, following the relations of Eq. (|3(J[1 . the numbers j\fj l of the corresponding operator 0jy . are such 
that, 

AL C = 2AS c = -AN + 2L c ^ 1/2 = ^ (QA^; 

AL S = 2AS s = A(N^-N l ) + 2L s> _ 1/2 = - J2 ^Q^Lh- (31) 

lc,l,=±l 

The first relation of Eq. (|20|l just states that the difference in the number of electronic creation and annihilation 
operators of the original A/"-electron operator O 1 ^ j equals the value AN of the electron number deviation generated 
by such operator. Similarly, the first relation of Eq. IjHlJI states that the difference in the number of rotated-electron 
annihilation and creation operators of 0jy ■ equals twice the value AS C = At] of the 7/-spin value deviation generated 
by that A/"-rotated-electron operator. Similar considerations apply to the second relations of Eqs. 112011 and 

We emphasize that all matrix elements of the general spectral-function expression (|29|l refer to regular energy 
eigenstates. Indeed, by changing from the spectral function representation J7J to (|29|l we have eliminated the explicit 
presence of —1/2 Yang holons and —1/2 HL spinons. This was done by evaluation of the contribution of these quantum 
objects to the A/"-electron spectral weight. Such a procedure corresponds to the computation of the commutator 
[n a =c. s (Sa) La -- 1/2 , O l Utj ] on the right-hand side of Eq. 

Our next step is the expression of the operator 0^ ^ for the general spectral-function expression l|29[) in terms 
of rotated-electron creation and annihilation operators. Here we use the results of Refs. 0, 0, 0] concerning the 
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expression of the rotated electrons in terms of ±1/2 holons, ±1/2 spinons, and cO pseudofermions. It is this direct 
relation that makes convenient the rotated-electron expression for the A/"-electron spectral functions. The expression 
of the local TV-electron operator Gjy j in terms of rotated-electron creation and annihilation operators is obtained by 
use of the following relation, 

oo 

< : >A-. ; < S (->a-..,' S «V.., ' X\ ^.V.,: 3 = l,2,-,N a ; l = ±l. (32) 

i=l 

Here S is the operator defined by Eqs. (21)-(23) of Ref. 0] and Q l jy j has the same expression in terms of rotated- 
electron creation and annihilation operators as 8^- ^ in terms of electronic creation and annihilation operators and 
thus Mo = M. It is given by, 

® l ^ 3 =V\U/t)e^ tJ V(U/t) = e- S 6^^e § ; j = 1,2, N a ; I = ±1 . 

The operators Oj^ • on the right-hand side of Eq. I|32l) such that i — 1,2, ... can be written as a product of Mi 
rotated-electron creation and annihilation operators and the value of the coefficient c\ is a function of n, m, and U/t 
such that c\ — > as U/t — > oo. For instance, for i = 1 and i — 2 we find, 

V^ek.j^^.QVo,,-]; j = l,2,...,N a ; l = ±l, (33) 

and 

^2®H 2 , 3 = \[SASi® l ^ 3 ]]\ j = l,2,...,N a ; l = ±l, (34) 

respectively, and the i > 2 operator terms are easily generated and involve similar commutators. For simplicity, here 
we omit the longer expressions of the latter terms. 

It is useful for the study of the spectral-function expressions to divide each CPHS ensemble subspace in a set of 
well-defined subspaces. The number deviation AN al , for the av = cO, si branches and the number N av = AN av for 
the av ^ cO, si branches can be expressed in terms of other related numbers as follows, 



AN au = AN F U + AN»f ; AN F V = AN^ +1 + AiV^ ; 2AJ F V = AN F V , +1 - AJ\£, _ x ; 
AiV^ L = AN°l F L + t Q° a j2n ; AN° a f = AN° a F +1 + AN° a > v F -i = AJ\£ ; av = cO, si o = ±1 ; 

N av = N F V + NZ F 5 N F V = N F ^ +1 + N F V< _! ; 2 j£ = +1 - <^ _ x ; av £ cO, si , (35) 

Here AN F V ±1 is the deviation in the number of = cO, si pseudofermions at the right (+1) and left right (— 1) 
Fermi points, AN F V and AJ F U are the corresponding number and current number deviations, respectively, and AN^f 
gives the deviation in the number of av = cO, si pseudofermions away from these points. Moreover, AN a ' u ±1 is the 
actual number of av pseudofermions created or annihilated at the right (+1) and left right (—1) Fermi points and 
Qq„/2 is the scattering-less phase shift on the right-hand side of Eq. Qllfl. For the av ^ cO, si branches, N F V t 
is the number of av pseudofermions with limiting bare-momentum value q = Lq^ v such that u = ±1, J F V is the 
corresponding current number, and N/// F is the number of av pseudofermions whose bare-momentum values obey the 
inequality \q\ < q° al/ . 

Let us also consider the number Np'/; nf of finite-momentum and finite-energy av pseudofermion particle-hole 
processes (A), which refers to the av = cO, si branches only 0. N^ NF is zero or a positive integer such that 
jyphNF _ \j\f av _ |AA^ F |]/2. Here M av gives the number of av pseudofermion creation and annihilation operators 
involved in the expression of the generators of the elementary processes (A). 

The J-CPHS ensemble subspaces are the subspaces of a CPHS ensemble subspace spanned by the excited energy 
eigenstates with the same values for the numbers Nj£ NF , N^ NF , AN F +1 , AN F 0> _ l5 AN F 1>+1 , AN F ^, and sets 
of numbers {N F V +1 } and {N F V _ x } for the av ^ cO, si branches with finite pseudofermion occupancy in the CPHS 
ensemble subspace. 

Use of Eq. H32fl for j = in the general spectral- function expression (|29|l leads to, 

oo 

B l u (k, uj) =J^4 E | E B l >\k,u)]; c'=l, l = ±l, (36) 

i=0 {AAT„, y },{L Qi _ 1/2 } {N>£ NF },{ANg„ ,},{Ng u J 
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where the summations J2{AN al ,} {l 1/2 } anc ^ J2{n pIiNF } {an f } {n f } run over CPHS ensemble subspaces and 
the corresponding J-CPHS ensemble subspaces of each of these spaces, respectively. The function B l,l (k, u>) on the 
right-hand side of Eq. (|3l))) reads, 

sW(fc, W ) = (J] ±-) J2 N a\(f-L; JC\Q l ^ \GS)\ 2 

a=c, s a f 

x d(uj~lAE^ cph ^5 k j Akj _ cphs ; l = ±l; i = 0,1,2,..., (37) 

where the summation runs over the excited energy eigenstates \f.L; JC) which span each J-CPHS ensemble 

subspace. Thus, there is a function B l ' l (k, uo) for each of these subspaces. (We recall that \f.L; JC) is the LWS of a 
state \f; JC) related to it by the general equation iffi^ .l 

Finally, let us use a notation for the number of spin-down and spin-up rotated-electron creation and annihilation 
operators of the operator 0^-. , such that i = 0, 1, 2, ... similar to that associated with the numbers Sf\ c ls of Eq. JHOJl. 
Thus, we introduce the numbers, 

Mi= tf?Js> ' = ±li * = 0,1,2... ) (38) 

l c , l s =±l 

which refer to the operator 0^-. y Here A/"/' 1 ; is the number of rotated-electron creation and annihilation operators 
for l c = — 1 and l c = +1, respectively, and with spin down and spin up for l s = —1 and l s = +1, respectively. The 
operator 0jy o ■ of Eq. (|3"2*|l has the same four rotated-electron numbers {Nfj ; a } = {^i c j s } as the corresponding 
operator 0^- ■ in terms of electrons. 

IV. ROTATED-ELECTRON SUM RULES, SELECTION RULES, AND ELEMENTARY PROCESSES IN 

TERMS OF PSEUDOFERMION OPERATORS 

In this section we provide sum rules and selection rules which for the PS arise from the direct relation between 
rotated electrons and the holons, spinons, and pseudofcrmions. Furthermore, we use such a relation to express the 
elementary rotated-electron processes in terms of the pseudofermion creation and annihilation operators. 

An important symmetry is that all six generators of the 77-spin and spin SU(2) algebras are invariant under the 
electron - rotated-electron unitary transformation 8\. Thus, the number of spin-projection a electrons equals the 
number of spin-projection a rotated electrons. This also applies to the deviations AS C and AS S in the 77-spin and 
spin values, respectively, generated by application onto the ground state of a A/"-electron operator. This symmetry 
implies that all i = 0, 1,2, ... operators 0^-. ^ on the right-hand side of Eq. Q32JI generate the same deviations AS C 

and AS S , as the operator 0^- ■ on the left-hand side of the same equation. It follows that the values Af { ' z ; for all 
these i = 0, 1, 2, ... operators with the same value of I = ±1 obey the following sum rules, 

AL C = 2AS c = -AiV + 2L c ,_ 1/2 = ^ (Qtffj, 5 

AL S = 2AS S = A{N ] -N i ) + 2L s ,_ 1/2 = - £ {l c l s )JV^ h ; l = ±l; i = 0,l,2,.... (39) 

; c ,z 3 =±i 

These rules provide useful information about the expression of all i — 1,2, ... operators 0^-. j. In addition to the 

rotated-electron creation and annihilation operators of 0jv" o ,ji such an expression includes pairs of rotated-electron 
creation and annihilation operators with the same spin projection a. Thus, such additional creation and annihilation 
operators only generate rotated-electron particle-hole excitations and do not change the net number of spin-projection 
a rotated electrons created or annihilated by application of the operators 0^. j onto the ground state. The general 
situation refers to A/"-electron operators that are not invariant under the electron - rotated-electron unitary transfor- 
mation. (The problem is trivial for those that are invariant, once 0^. , = for i > in that case.) The precise form 
of 0^-. j for i = 1,2, ... depends on the specific A/"-electron operator under consideration. However, a general property 
that follows from the relations 139|) is that for increasing values of i = 1, 2, ... the operators 0^. j are constructed by 
adding to 0^ o • an increasing number of particle-hole elementary spin-projection a rotated-electron pairs. 

In equation (|32|l the A/"-electron operator 0^ j is expressed in terms of rotated-electron creation and annihilation 
operators. From the direct relation between the rotated electrons and the holons, spinons, and pseudofcrmions it is 
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straightforward to find useful selection rules. Such rules refer to restrictions in the values of the number of —1/2 
holons and thus of 2^-holon composite cv pseudofermions generated by application onto the ground state of each of 
the i = 0, 1, 2, ... operators ®\f it j of expression [The expression of these operators determines the value of the 

TV-electron spectral function of Eq. (|36() . as confirmed by the form of the related functions l|37|) .] For the PS excited 
energy eigenstates which have finite overlap with the TV-electron excitations, the values of the — 1/2-holon number, 
M c ,-1/2, and number of finite-energy and finite-momentum cO pseudofermion particle-hole processes, N^q NF , of the 
elementary processes (A) Q are restricted to the following ranges, 

oo 

u=l l s =±l 

N% NF = 0,1,.., mini ]T tf l -i t i s , E i i = 0,l,2,.... (40) 



Here the numbers A/"/'*; are those of Eq. I|38(l for the operator Oj^ y The i = operator 0jy o j has the same 
expression in terms of rotated-electron creation and annihilation operators as the corresponding operator 0^- • of Eq. 
(|32|l in terms of electronic creation and annihilation operators. Therefore, for i = the selection rules given in Eq. 
|@DJ read, 

^c,-i/2=i c ,-i/2+£^A ct/ = 0,l,..., tf-i,h-> N^ NF =0,1,..., mini £ AfL 1Js , £ \ , (41) 

V=\ i„=±l U s =±l ( a =±l J 

where the numbers l = A7'l' < j ) l are those of Eq. (|30[1 specific to the operator Ojy .. 

The first exact ground-state charge selection rule of Eq. (|40|l concerning the number of —1/2 holons, M c , _i/2, is 
equivalent to the following selection rule involving the number deviation — AM C — AM, = A7A c o, 

E < AM C < Mi ; —Afi < AA c0 = AM S < - E ' * = °> 2 ' - ■ ( 42 ) 

Indeed, the combination of the inequalities (|42|) with the relations (27)- (29) of Ref. [8j and that of Eq. (|39|l . readily 
confirms the equivalence of the first selection rule given in Eq. (|40|l and that of Eq. I|42|l . 

Moreover, the —1/2 spinon number deviation AM s _ 1 / 2 is fully determined by the value of the —1/2 holon number 
of the first selection rule of Eq. l|40|l and reads, 

oo 

AM s> _ 1/2 = L a ,_ 1/a +AN sl + J2» N ™ = AN l ~ M c,-i/2 • (43) 

Equations (THJ, and ® define sum rules obeyed by the values of the deviations in the quantum-object 

numbers and Eqs. I|21|) . (|40|) . and (|42|l correspond to selection rules for the permitted values of these deviations. Such 
sum rules and selection rules define the set of CPHS ensemble subspaces which contain the excited energy eigenstates 
with finite overlap with the A/"-electron excitations under consideration. 

While the above rules are exact, direct evaluation of the weights by the method introduced in Sec. V and further 
developed in Ref. Q reveals that 94% to 98% of the A/"-electron weight corresponds to excited energy eigenstates with 
numbers in the following range, 



oc 



L^-i/a+E^- 1 )^^ ' 1 ' 2 '-' E 5 ^is^ l it,i s \ * = 0,1,2,.... (44) 
For i = the relation l|44fl can be written as, 

oo 

L.,_i/ a + X>-l)iV», = 0,l,2,..., s i«i.K,i.> (45) 

where the numbers 7V/ c ls — tyfi are those of Eq. (|3*U|l specific for the operator 9^ • of Eq. ijS^ . 

Local —1/2 holons (and —1/2 spinons) correspond to local 2^-holon composite cv pseudofermions (and 2z/-spinon 
composite sv pseudofermions). Local av pseudofermions are associated with the operators /1 jjQJ/ and f Xj ,au on the 
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right-hand side of Eq. (34) of Ref. Let us denote the rotated-electron spin projections a =|, | by er = —1/2, +1/2, 
respectively, and consider the elementary processes of the A/"-electron excitations in terms of occupancy configurations 
of local ±1/2 holons, ±1/2 spinons, and cO pseudofermions [l7j : 

(i) To create one spin-projection a = ±1/2 rotated electron at the unoccupied site j, we need to annihilate a local 
±1/2 holon and create a local cO pseudofermion and a local ±1/2 spinon at the same site. Annihilation of a spin- 
projection a = ±1/2 rotated electron at a spin-projection a = ±1/2 rotated-electron singly occupied site j, involves 
the opposite processes. 

(ii) To create one spin-projection a — ±1/2 rotated electron at a spin-projection a = ±1/2 rotated-electron singly 
occupied site j, we need to annihilate a local +1/2 spinon and a local cO pseudofermion and to create a local —1/2 
holon at such a site. Again, to annihilate a spin-projection a = ±1/2 rotated electron at a rotated-electron doubly 
occupied site j, involves the opposite processes. 

(iii) The creation of two rotated electrons of opposite spin projection onto the unoccupied site j involves the 
annihilation of a local ±1/2 holon and the creation of a local —1/2 holon at such a site. Annihilation of two rotated 
electrons of opposite spin projection onto the doubly-occupied site j, involves the opposite processes. 

(iv) The annihilation of one spin-projection a — ±1/2 rotated electron and creation of one spin-projection a = ±1/2 
rotated electron at the singly-occupied site j, involves the annihilation of one local ±1/2 spinon and the creation of 
one local +1/2 spinon. 

Other processes can be expressed as suitable combinations of the above elementary processes. The local rotated- 
electron operator terms which transfer spectral weight from the ground state to each of the J-CPHS ensemble subspaces 
appearing in the state summation of the spectral-function expression (|3*S|) have a specific and uniquely defined form in 
terms of cO pseudofermion and composite ctv pseudofermion creation and annihilation operators. In order to find the 
pseudofermion form of these operator terms it is crucial to take into account the initial ground-state pseudofermion 
occupancies, given in Eqs. (C.24)-(C.25) of Ref. [8|. (We recall that the pseudoparticle-number values of the latter 
equations equal those of the corresponding pseudofermion numbers.) 

Before illustrating how the elementary processes (i)-(iv) are generated by the pseudofermion creation and annihi- 
lation operators, it is convenient to provide some basic rules for the use of the latter operators. Since following the 
use of the relations of Eq. all matrix elements are between the ground state and regular excited states, in the 
processes considered below, the deviations in the numbers of Yang holons (a — c) and HL spinons (a — s) are such 
that AL a = AL a+ i/ 2 - Some of these processes involve creation or annihilation of +1/2 Yang holons and/or +1/2 
HL spinons. However, we recall that within the pseudofermion representation, the +1/2 Yang holons and +1/2 HL 
spinons do not appear explicitly. Such processes are taken into account by the deviations in the number of discrete 
bare- momentum (and canonical-momentum) values and effective av lattice sites of the cv ^ cO branches and/or sv 
branches, respectively |17| . Given the values of the corresponding pseudofermion number deviations, this is readily 
confirmed if one compares the number (B.6) of Ref. Q of discrete bare-momentum values and of effective lattice 
sites of the excited energy eigenstate and ground state CPHS ensemble subspaces. Since AL a — AL a+1 / 2 , note 
that following Eqs. (B.6) and (B.7) of Ref. 0, the value of the number N*„ changes when the value of the number 
L c = L c , -1-1/2 of +1/2 Yang holons and/or L s = L s+1 / 2 of +1/2 HL spinons also changes. Thus, creation and 
annihilation of +1/2 Yang holons (and +1/2 HL spinons) are processes that are taken into account in the definition 
of the effective cv ^ cO pseudofermion lattices of the initial ground state and excited energy eigenstates. 

In the following we provide different examples of local rotated-electron operator expressions in terms of pseud- 
ofermion creation and annihilation operators. For simplicity, each of such pseudofermion expressions corresponds to 
the term of the local rotated-electron operator which transfers spectral weight from the initial ground state onto a 
single excitation J-CPHS ensemble subspace. Such a pseudofermion term includes a coefficient factor 1/Cj whose 
value is well defined for each subspace. The full local rotated-electron operator term which transfers spectral weight 
from the ground state to a given J-CPHS ensemble subspace is the product of that studied here by another pseud- 
ofermion operator term given in the ensuing section. (The general form of 1/Cj is also given in that section.) The 
latter operator transfers from the ground state to the J-CPHS ensemble subspace the part of the spectral weight 
which corresponds to the processes (B) and (C), whereas the pseudofermion terms studied here transfer the part of 
that weight associated with the processes (A) . 

A i — local rotated-electron operator Ojy o • always has one or a few dominant CPHS ensemble subspaces which 
correspond to the whole spectral weight transferred from the ground state in well defined limits. For one-electron 
problems such that A/"o = Af — 1 this refers to the limits where the spectral-weight distribution is 5-function like, 
as in Eq. (77) of Ref. For Afo = Af = 2, to the limits where such a distribution can be expressed as a simple 
integral whose integrand is a 5 function, as in Eq. (78) of the same reference. In the general Afo = Af > 2 case, to the 
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limits where the spectral-weight distribution can be written as an integral whose integrand is a product of [TV — 1] 
5 functions. For instance, for TV = 1 this occurs for U/t — ► 0. For the one-electron problem the amount of spectral 
weight transferred from the ground state to the set of J-CPHS ensemble subspaces contained in the dominant CPHS 
ensemble subspaces is weakly dependent on U/t: while for U/t << 1 it corresponds to the whole spectral weight 
transferred from the ground state by the local rotated-electron operator, for U/t >> 1 it corresponds typically to 
more than 0.94% of that weight. (Comparison of the amount of transferred weight for U/t << 1 with that obtained 
by use of the methods of Ref. [ljj for U/t >> 1 confirms such a weak U/t dependence.) For U/t finite there arise 
an infinite number of pseudofcrmion terms, each corresponding to a J-CPHS ensemble subspace compatible with the 
local rotated-electron operator. Another example is the TV = 2 charge dynamical structure factor, where there are 
different dominant CPHS ensemble subspaces for U/t — > and U/t — > oo, respectively. In this case the amount of 
spectral weight transferred from the ground state to the set of J-CPHS ensemble subspaces contained in the dominant 
CPHS ensemble subspaces is a decreasing (and increasing) function oiU/t for the U/t — > (and U/t — > oo) dominant 
subspaces and vanishes as U/t — > oo (and U/t — > 0). Again, for intermediate finite values of U/t there arise an 
infinite number of pseudofermion terms, each corresponding to a J-CPHS ensemble subspace compatible with the 
local rotated-electron operator. However, for alH = local rotated-electron operators 0^ o ^ the pseudofermion terms 
associated with the dominant subspaces together with a small number of other terms correspond to more than 99% 
of the spectral weight. It follows that in applications of the pseudofermion dynamical theory introduced here and in 
Ref. 0] only a finite number of pseudofermion terms should be considered. 

The av — cO, si pseudofcrmion number deviations and av ^ cO, si pseudofcrmion numbers are related to the 
rotated-electron number deviations by Eqs. I|18|) and (|19|) . Given the values of the av = cO, si pseudofermion number 
deviations and av ^ cO, si pseudofermion numbers of the specific J-CPHS ensemble subspace under consideration, 
the expression of the local rotated-electron operator in terms of pseudofermion creation and annihilation operators 
is always uniquely defined. Let us start by providing some of the simplest pseudofermion operator terms of local 
rotated-electron operators. For local one- and two-rotated-electron operators these operator terms involve in general 
av pseudofermion creation and/or annihilation operators belonging to branches such that v < 2. The case of other 
terms associated with excitation J-CPHS ensemble subspaces generated from the ground state by processes involving 
creation of composite av pseudofermions for v > 1 is discussed later. 

In the following expressions the av effective lattice integer site index j' is such that j' — 1, 2, N* v 17]. An 
important property is that an operator whose expression in terms of rotated-electron operators is local at Xj = ja = j 
can be written as a product of local av pseudofermion operators at Xy ~ Xj where Xj> = j' or av . Here and in all 
expressions given below j' is defined for the av =^ cO branches as the closest integer number to jn* ai/ , whereas j 1 = j 
for av — cO. We note that for the former branches the site j' occupied by one av pseudofcrmion corresponds to 2v 
sites of the rotated-electron lattice. Thus, \xy — Xj\ is always smaller than the very small intrinsic uncertainty which 
corresponds to the 2v rotated-electron lattice sites occupied by the local av pseudofcrmion. Moreover, we emphasize 
that the rotated-electron lattice site j associated with the effective av ^ cO lattice site f « jn* av defined above always 
belongs to the domain of 2v rotated-electron lattice sites of f . Here and below we use the equality j' — jn* alJ to 
denote the integer number j' defined above. Thus, the site j 1 is such that j 1 = j for av — cO operators, j 1 = jn-\ for 
av = si operators, j' = j[l — n] for cv ^ cO operators when n < 1, and j' = j[n^ — for sv ^ si operators when 
to > 0. The following local rotated-electron operator expressions in terms of pseudofermions, whose coefficient Cj is 
different for each operator, refer to the elementary processes (A) subspace: 

(i) One of the simplest processes for creation of one spin-down rotated electron at the unoccupied site j involves 
the creation of a local cO pseudofermion with the operator /J. c0 and of a local si pseudofermion with the operator 

fx jlt sl such that f =J n h 

4,i( 1 -^t) = ^:4„ s i4,co. (46) 

and thus Ej t i (1 — f) = fx^cO fxj,,si refers to annihilation of one spin-down rotated electron at the singly- 
occupied site j. In Eq. (I4(j|) . hj >a — Cj a Cj. cr is the local spin-projection a rotated-electron density operator. 

To create one spin- up rotated electron at the empty site j, a simple process corresponds to create a local cO 
pseudofermion with the operator /J. c0 , 

4,T( 1 -^i) = ^4,cO> (47) 

and Cj j i (1 — fij^i) = -jk- fx^cO to annihilation of one spin-up rotated electron at the singly-occupied site j. Such 
processes also involve creation and annihilation, respectively, of an empty site in the effective si lattice. When 
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the initial ground state belongs to a m = CPHS ensemble subspace, there is for the former process a single si 
pseudofermion hole in the excited state, which corresponds to the created site. 

We note that the processes of the first expressions of Eqs. (|46|l and l|47(l also involve the annihilation of a +1/2 
Yang holon, whereas the processes of the second expressions of the same equations involve the creation of a +1/2 Yang 
holon. Similarly, the processes of the first and second expressions of Eq. (|47|l involve the creation and annihilation, 
respectively, of a +1/2 HL spinon. In the remaining cases considered below we do not specify the elementary processes 
of creation or annihilation of +1/2 Yang holons and +1/2 HL spinons, which are taken into account implicitly by the 
pseudofermion description, as discussed above. 

(ii) One of the simplest processes associated with the creation of one spin-up rotated electron at a spin-down 
rotated-electron singly occupied site j involves the annihilation a local cO pseudofermion with the operator f x , c o and 
of a local si pseudofermion with the operator f Xj ,,si and the creation of a local cl pseudofermion with the operator 

fl „,ci sucn that f — 3 n \ and j" = j[l — n], respectively, 



c). T fij, i = ^- ft.„ t cl f Xj , co f Xj , ,si- (48) 

Then £j. j hj t j = (jjft., s i f x . c o fx } n,ci refers to annihilation of one spin-up rotated electron at a doubly occupied 
site j. Moreover, to create one spin-down rotated electron at a spin- up rotated-electron singly occupied site j, a 
simple process corresponds to annihilate a local cO pseudofermion with the operator f x , c o and to create a local cl 
pseudofermion with the operator / J , cl such that j' = j[l 



cl,4»*i,t = cj4,,ci./W,cO. (49) 



In this case Cj^hj^ — -^j-f Xj co/a^.ci corresponds to annihilation of one spin-down rotated electron at a doubly 
occupied site j. 

(iii) A simple process involved in the creation of two rotated electrons of opposite spin projection onto the empty 
site j corresponds to creation of a local cl pseudofermion with the operator f^ r cl such that j' = j[l — n], 

4l<T = ^4,ci- (50) 

It follows that Cj,]Cj,i — fj- f Xj ,,ci refers to annihilation of two rotated electrons of opposite spin projection onto 
a doubly-occupied site j. This involves annihilation of a local cl pseudofermion with the operator f x ., tC i such that 
f = j[l - n}. 

(iv) One of the simplest processes associated with the annihilation of one spin-up rotated electron and creation of 
one spin-down rotated electron at the singly-occupied site j, involves the creation of a local si pseudofermion with 
the operator /J ; sl such that j' = jn-f, 

£},i^T = ^4,»i- ( 51 ) 

Then cj j Cj. T — cT7 corresponds to annihilation of one spin-down rotated electron and creation of one spin-up 

rotated electron at the singly-occupied site j. This involves the annihilation of a local sl pseudofermion with the 
operator f x ., lS i such that j' = jn\. 

For local rotated-electron operators generating more complex processes involving creation or annihilation of sev- 
eral rotated electrons, the creation and annihilation of local cO pseudofcrmions is always associated with creation 
and annihilation of rotated-electron singly occupied sites, respectively. Since the local cO pseudofermions and cO 
pseudofermion holes occupy the same sites ji and jh as the rotated-electron singly occupied sites and rotated-electron 
doubly-occupied and unoccupied sites, respectively, there is a one-to-one correspondence between the rotated-electron 
and cO pseudofermion algebras. Here we have used the site notation of Ref. |l7j |. 

However, once the composite local av ^ cO pseudofermions have internal structure that involves 2v rotated-electron 
sites with different index j, the operational relation of rotated electrons to such composite quantum objects is more 
involved 01- This justifies why the expressions of the local rotated-electron operators in terms of creation and 
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annihilation pscudofcrmion operators involve a superposition of different pseudofermion expressions, corresponding to 
the set of compatible J-CPHS ensemble subspaces. Nevertheless, creation onto the ground state of a cv pseudofermion 
(and sv pseudofermion) always involves v rotated-electron doubly occupied sites and v unoccupied sites (and v 
spin-down rotated-electron singly occupied sites v spin- up rotated-electron singly occupied sites). In general, the 
excited-energy-eigenstate v rotated-electron unoccupied sites (and v spin-up rotated-electron singly occupied sites) 
of a created cv pseudofermion (and sv pseudofermion) are generated by annihilating an equal number of +1/2 Yang 
holons (and +1/2 HL spinons) of the initial ground state. 

For the creation of a local cv pseudofermion, each of the v new created rotated-electron doubly occupied sites can 
result from creation of a rotated-electron pair onto an unoccupied site or of a rotated electron onto a singly-occupied 
site. The latter case involves always one of the elementary processes associated with the pseudofermion terms given 
just after Eqs. (|46|l and (|47|l . On the other hand, for the creation of a local sv pseudofermion such that v > 1, the v 
involved spin-down rotated-electron singly occupied sites can result from creation of spin-down rotated electrons onto 
unoccupied sites or from recombination of pre-existing ground-state si pseudofermions, as further discussed below. 

For instance, let us consider two J-CPHS ensemble subspaces contained in different CPHS ensemble subspaces 
which except for the occupancies of the cl and c2 branches have the same pseudofermion numbers. For such branches, 
one has {N c \ — 2, AT c2 = 0} for the J-CPHS ensemble subspace (I) and {N c i = 0, N C 2 = 1} for the J-CPHS ensemble 
subspace (II) . Let us consider that the local rotated-electron operator behind the transitions to both subspaces is the 
same and involves creation of two rotated electrons of spin projections a =| and a =1 onto the spin-down singly 
occupied site j and spin-up singly occupied site j + 1, respectively. For the subspace (I), this corresponds simply to 
the process of Eq. I|48|l for the site j and the process of Eq. I|49|l for the site j + 1. For the subspace (II), in order 
to create two rotated electrons of spin projections a —\ and a —I onto the spin-down singly occupied site j and 
spin- up singly occupied site j + 1, respectively, we need to annihilate two local cO pseudofermions with the operators 
fxj,cO and f Xj+1 ,cO and a local si pseudofermion with the operator f Xj ,,si such that j 1 = j'nj and to create a local c2 
pseudofermion with the operator /|.„ c2 sucn that j" = iW ~ n L 

5 iT 5 Il^i^T = -^■ ] fl ] „,c2fx ] ,cofx 3 ,,sl fx J + 1 ,cO- (52) 

It should be mentioned that in spite of the annihilation of one si pseudofermion, this process does not involve the 
corresponding creation of a si pseudofermion hole. Indeed, it involves the annihilation of the site j' = jrij in the 
effective si lattice. Thus, when the initial ground state belongs to a m = CPHS ensemble subspace, in spite of the 
annihilation of the si pseudofermion the excited state corresponds to a fully occupied si band, as the initial ground 
state. 

A similar process gives rise to creation of a local s2 pseudofermion provided that creation of the two rotated-electron 
doubly-occupied sites is replaced by creation of two spin-down rotated electron singly occupied sites. However, in this 
case there is the possibility that one (or both) the spin-down spinons needed for creation of the local s2 pseudofermion 
is (or are) generated from annihilation of one (or two) ground-state si pseudofermion(s). Such processes can dress 
any rotated-electron process and are behind the occurrence of an infinite number of compatible J-CPHS ensemble 
subspaces for each local rotated-electron operator. These non-dominant pseudofermion processes do not obey the 
relation (|45|l (which is not an exact rotated-electron selection rule) and for all finite values of U/t amount to less 
than 6% of the rotated-electron spectral weight. For instance, in order to create one spin-down rotated electron at 
the empty site j, in addition to the pseudofermion process (|46|l there is for instance a process corresponding to the 
creation of a cO pseudofermion with the operator c0 and of a s2 pseudofermion with the operator fl u s2 such 

that j" = j[n^ — njj and to the annihilation of a si pseudofermion with the operator fl r sl such that j' = jn^, 

sJa f 1 - fi i,r) = ^; fl r ^2 /«y,.i 4, c0 • (53) 

We emphasize that the amount of spectral weight transferred from the ground state by the operator (|46|l is much larger 
than that transferred by the operator Indeed, the J-CPHS ensemble subspace associated with the expression l|46|) 
belongs the dominant CPHS ensemble subspace of the local rotated-electron operator, whereas the J-CPHS ensemble 
subspace corresponding the expression l|53|) does not. Note that in the present case the two competing J-CPHS 
ensemble subspaces can have the same pseudofermion deviation numbers and values except for the occupancies of the 
sl and s2 branches. Thus, the processes generated by the operators (|46H and (|53|) correspond to J-CPHS ensemble 
subspaces belonging to different CPHS ensemble subspaces such that {AiV s i = 1, N S 2 = 0} and {A7V s i = —1, iV S 2 = 
2}, respectively. Similar dressing processes involving creation of sv pseudofermions belonging to v > 1 branches by 
annihilation of one to v ground-state sl pseudofermions can occur for all rotated-electron processes but correspond 
to very fast decreasing values of the amount of spectral weight transferred from the ground state for increasing 
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number of pseudofcrmion processes. Moreover, we recall that the subspace summation on the right-hand side of the 
spectral-function expression l|36ll is limited to the compatible CPHS ensemble subspaces: their pseudofermion number 
deviations and numbers obey the sum rules ifT ^I . ffi ty . and and selection rules lf2"T |) . pi jl . and 

We could present here other pseudofermion terms of increasing complexity, corresponding to the local rotated- 
electron operators considered above. However, the amount of spectral weight transferred from the ground state 
by the pseudofermion operator terms describing the above-mentioned dressing processes involving creation of sv 
pseudofermions such that v > 1 by annihilation of an increasing number of ground-state si pseudofermions decreases 
very rapidly for increasing values of v. Also the spectral weight transferred from the ground state by the operator 
terms with increasing value for the index i of the expression 1M2I) of the general operator • decreases very rapidly. 
For instance, for the one-electron spectral weight the contributions from dressing sv pseudofermion processes for 
v > 2 and the terms of index i > 1 of the expression (|32fl for M = 1 are typically beyond numerical measurability. 
Therefore, as far as numerical measurability is concerned, only a few pseudofermion terms contribute to the actual 
electronic spectral- weight distributions 0. 



V. THE SPECTRAL FUNCTION AS A CONVOLUTION OF PSEUDOFERMION SPECTRAL 

FUNCTIONS 

In this section we express the spectral functions <|37[1 as a convolution of pseudofermion, independent —1/2 holon, 
and independent —1/2 spinon spectral functions. The excited energy eigenstates appearing on the right-hand side of 
Eq. (|r?T|> can be written as the following pseudofermion Slater determinant, 

\f-L,C)=X\F} av \Q)- F} av = \{ [Afavi&flw + ll-tfauiqj)]]. (54) 

Here and in other expressions below |0) is the pseudofermion vacuum such that fq- a v\fy = for all av branches and 
canonical-momentum values. In turn, according to Eqs. (C.24) and (C.25) of Ref. |8j], the ground state corresponds to 
a canonical- momentum densely packed occupancy for the cO and si pseudofermion bands and the Slater determinant 
has the following simpler form, 



aU=cO, Sl ^j—~1°Fav 

qFav,+l 9Pou, +1 

F^-qs. olv = II 4-5 F$_ GSiav = n /Ui <w = cQ,*l. (55) 



The generators F^_ GS av and Fj_ GS av given here correspond to densely packed distributions introduced below and the 
discrete canonical-momentum values of the pseudofermion operators f\^ av of their expressions are those of the CPHS 
ensemble subspace which the ground state \GS) and the excited state \f.L; C) of Eq. l(5~4*|) belong to, respectively. 
The Fermi points appearing in the products of their expressions of Eq. read qFuv, ±1 = ^QFau ^ fa / L]AN^ ±1 
and qFav, ±i = ^1% av ± fa / L][AN£ U> ±1 ±Q® lf (±q Fai/ ) /2n], respectively, where the deviation numbers AJV°^f ±1 and 
AiV ( f !/) ±1 are those of Eq. ffify . 

The excited-energy-eigenstate canonical- momentum distribution function N av {q~j) on the right-hand side of Eq. 
(|54|l can be written as, 

Nav {q~j ) = + AAC/fe) ; A^fe-) = A^fe) + AAJ^fe) I <*" = cO, sl 

ACfe) = ^feO+AACfe); ai/ = cO, sl; M av (q j )=AAf a l F (q j ) + ^ u (q j ); av cO, sl . (56) 

Here N~„(qj) and N~®(qj), such that M~y{qj) — N~®(qj), correspond to the excited densely packed distributions 
I~Iqi/— cO, sl -GS,aJ0) an< ^ riQi/=cO,si Fl GS a jO), respectively. Furthermore, the ground-state distribution J^auilj) 
is both such that q~j = qj and A/2 1/ (g J ) = NjL,(qj), where N® u (qj) is the ground-state bare-momentum distribution 



function given in Eqs. (C.1)-(C.3) of Ref. \$. Thus, AA/"£,fe) = AN£ v (qj) describes av = cO, sl pseudofermion 
addition to or removal from the Fermi points. Moreover, AJ\f^ u (<lj) — AN^f (qj) describes av pseudofermion creation 
and/or annihilation away from the Fermi points for the av = cO, sl branches and creation of av pseudofermions at 
canonical-momentum values such that \q~j\ < q® av for the av ^ cO, sl branches, whereas AAf^(qj) — AN^(qj) 
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describes creation of av ^ cO, si pseudofermions at the limiting canonical-momentum values qj = ±q av - Finally, 
the deviation AAfP^ F (qj) — AN^ F (qj) corresponds to low-energy and small-momentum av = cO, si pseudofermion 
particle-hole processes. For n = 1 (and/or m — 0) the excitation subspace is such that AJ\f^ F (qj) — AN^ F (qj) = 
(and/or A7vj/ lF (^) = AN^ F (q 3 ) = 0). The above deviations are such that, 

+qL +iL +iL 

£ AN£ F %)=0; J2 Wf» F (q j ) = AN» F ; £ A^ F v (qj) = AN F V ; av = cO, si . 

9j=-1a U 9J=-9°„ 3i=-9°t. 

Since for the <w 7^ cO, si pseudofermion branches there is no occupancy in the initial ground state, the canonical- 
momentum distribution function is such that, =I g o AJ\f av (q~j) = N^ F + N^ v . [See Eq. l|35l) .] 
The generator f\ av given in Eq. I|54|) can be written as, 



/, cuv 

Ff t av = ^p-h, aw ^J-NF, av ^J-GS, av > ai/ = s l > 

= 4, tt ,4, ,; «^ c o, si. 

The expressions of the generators Fj_ NF av , F^_ h av , and F^ F av in terms of pseudofermion creation and annihilation 
operators are given in Eqs. (B.l), (B.4), and (B.5) of Ref. respectively, and that of the generator Fj_ GS av is 
provided in Eq. (|55|) . whereas the generator F F av creates av ^ cO, si pseudofermions at q = iq^ u - The generators 
^p-h ay Fj-GS av an d Fj-nf av are associated with the deviations Nav(Qj) and AN^/ F (qj) of Eq. JSflJ), respectively, 
Fp_ h av corresponds to the deviation AAf^ F (qj) of the same equation and thus generates the low-energy and small- 
momentum av — cO, si pseudofermion particle-hole processes (C). The operators Fjj F and F F 



, av 



av 7^ cO, si branches and are associated with the deviations AN/// v F (qj) an d ASf F v (<lj), respectively. 

The precise expression of the spectral functions of Eq. (|37[) depends on the specific form of the local operator 
Q\f. whose expression includes contributions from all av branches with finite pseudofermion occupancy in the 
corresponding J-CPHS ensemble subspace. For each such a subspace that operator expression and the corresponding 
coefficient Cj appearing in Eqs. |@SJ-(J53J have the following general form, 

©kj = ^-tn^a*]; C J = e ijAPj [Gc/Gj] ; Gj = JJ [Wa»)« v )^ + [1 - 0(M av )] 

OLV OiV 

H av = \AN» F \+2Ng rfF , au = cD,8l\ N av = N» F , av ? cO, si , (57) 

where the values of the positive real coefficient Gc and momentum deviation APj are given below and we recall that 
f = jn* au denotes the integer number closest to jn av . In equation (|57|l and in the remaining of this paper 9(x) is 
such that 9(x) = 1 for x > and 9(x) = for x < 0. 

The value of the coefficient Gc appearing in the general expression for the operator ©j^. j given in Eq. I|57|l is the 
same for all J-CPHS ensemble subspaces belonging to a given CPHS ensemble subspace, whereas that of Gj and thus 
of \Cj\ is specific to each J-CPHS ensemble subspace. Furthermore, when the expression of the local operator 0\f j 

of Eq. is independent of U/t, the same occurs for the related operators ©jy j and ©V , j 01 Eqs. (|27|l - lj2%)l and 
(|32|l , respectively. It follows that in the case of the i = operator 0^ ^ , the value of the coefficient Gc appearing in 
its expression given in Eq. (|57(l is also independent of U/t. Fortunately, for the dominant CPHS ensemble subspaces 
considered in the previous section, such a value can be found from analysis of the problem for U/t — or U/t = 00. 
One then finds that Gc = 1 in the ©V ,j — cfe'lllcw av] expression corresponding to the dominant CPHS 
ensemble subspaces. 

The operator Q l -J av appearing in Eq. (|57|l has the following general form for the av branches with finite pseud- 
ofermion occupancy in the J-CPHS subspace, 

oM . „-.AP°„j'«L al,NF,i a;,_f,,: 

^j'.av c ^j'tOtv ^av • 

Here the operators ©y^' and Q l ^ F ' 1 are associated with the elementary processes (A) and (B), respectively, and the 
phase-factor momentum lAP® y is given below. Considering that [xji — Xj] — 0, it is such that —i^2 av AP^,j'a'„ = 
—ijAPj. Here the summation Y] m , is over the av branches with pseudofermion occupancy in the J-CPHS ensemble 
subspace and APj is the momentum deviation of Eq. 1)57(1 . which reads, 

APj = 5>p°„. 
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Let us proceed by studying the phase factor e lAP aJ a '<^ and operators 0^'* and whose product gives the 



operator Qj? av . We start by characterizing for each av branch the elementary processes that originate the phase factor 
e -iAP al/ j a a „ _ ^p^g expression of the local operator 0^,' av does not involve the generator of the elementary processes 
(C), F£_ h av , but provides the momentum lAPg£ F for such processes through a phase factor, e~ %AP °^ F ^ a °". Here 

lAP'^f = 1[2tt/L] [m a v,+i — rn a v,-i] is the small momentum deviation generated by the av — cO, si pseudofermion 
particle-hole elementary processes (C) and m av> ±\ is the number of such processes of momentum ±[2ir/L] in the 
vicinity of the Fermi points ±q Fav . Furthermore, each av = cO, si pseudofermion created or annihilated at the Fermi 
point l q Fau by the elementary processes (B) contributes with a phase factor e~ u q F<*vi a av r e + " a ° „ ; respectively, 
where i = ±1. This leads to a phase factor e- i2q F<*» AJ %? j ' a °" , where 2AJ°f = AN% F +1 - AN^ F _ X and the number 
deviation AN®' F ±1 is that defined in Eq. i|35|) . Moreover, each cv ^ cO and sv ^ si FP scattering center created by the 
elementary processes (B) contributes with a phase factor e - " 9 ^-? a <=o and e~ uqp <' ^ Qrf) e +u2qFsi: > asi , respectively. On 
the other hand, the scattering-less bare- momentum shift contributes with a phase factor e~ l ^ a "' L ^ a ^ for each of the 
N® v pseudofermions of the initial ground state, what gives \e~ l ^ av ^ L ^ a <n>] N av = e~ %qFa "^^"l^ "™ with av = cO, si. 
Adding all these contributions leads to the above net phase factor e~ %APa »i a av f or the av — cO, si branches whose 
phase- factor momentum reads IAP® V = l[AP^ F + AP F V ]. For densities in the ranges < n < 1 and < m < n, 
the momentum AP F V appearing in that phase factor is given by AP F = 4k F [AJ F + Y^uLx ^cu + Y^=2 Jfv] an d 
APf x = 2kpi[AJ F 1 — 2jy^=2 Jsl] f° r the av = c ® an d av = s l branches, respectively. It results from the current 
contributions associated with the cO and si Fermi points ±2k F and ±k F [, respectively. Moreover, each cv =/= cO 
pseudofermion created by the elementary processes (A) contributes with a phase factor e^ l< - 1+v ^ 7 '^ a °" what leads to 
a net phase factor e~ lAP "^ a " v for the cv ^= cO branches such that AP C „ = (1 + v)itJ\[ C v Finally, the phase- factor 
momentum vanishes for the sv ^ si branches. Thus, the phase-factor momenta contributing to APj read, 

IAP° V = l[AP^ F + AP F V ] , av = cO, si ; IAP° V = 1(1 + v)nM c , , cv £ cO ; IAP° SU = , sv ^ si . 

Next we consider the operator Q l ^ Ft - For the cv ^ cO and sv ^ si branches, that operator creates 2vN F v 
independent —1/2 holons of momentum it and 2vNf v independent —1/2 spinons of momentum zero, respectively. (The 
only effect of the creation of the corresponding N F V cv ^ cO FP scattering centers and N F V sv ^ si FP scattering centers 
0, respectively, is the above contribution to the phase factor e~ lAP ° o: > a °° and e~ lAP ° o: > a "o e - zAP °ii a °i ; respectively.) 
For the av = cO, si branches the operator & l ^ Fl is such that, 

U{^ N U n [t 

l=±1 i'=l j"=l 

\ AN L J-l N' 



\ t =±i ,-'=i j"=i 

\ AN L J-i jv* 

+ 0(-Ai\O n [Ee*- 2117 ^-' 3 "^/,,,,.]}; ^ = cO,sl 



i'=Q j»=l 

Thus, that operator can be written as, 

AN L,, \ AN L,,\~ 1 



e 



n{e(AiV^J [] 4o„ +27ri7i) , a , + e(-AiV^J JJ / t( ^_ 27ri , /iW }; a^cO,sl 



t=±l " i'=l i'=0 



where we omitted corrections of order 1/L to the momentum value appearing in the phase factor. However, these 
corrections must be considered in the momentum of the pseudofermion operators. In the above two equations the 
pseudofermion operators correspond to spatial and canonical-momentum variables, respectively, and in the argument 
of the exponentials appearing in these equations and in other equations given below i is the usual imaginary number. 
(It is not the index i of the operator Oq^'\ whereas i' is a summation index.) Moreover, in the above equations and 
in the remaining of this paper O(cc) is such that Q(x) = 1 for x > and Q(x) = for x < 0. 

We finish the study of the operator O'', 4 av by considering the operator Q 1 ^ 1 ^ associated with the elementary 
processes (A), whose expression refers to a given J-CPHS subspace. In order to arrive to that expression in terms of 
the local av pseudofermion creation and annihilation operators, we recall that our study refers to spectral functions 
of form J7J whose operator O 1 / • expression involves Af elementary electronic operators which create or annihilate 
electrons in a compact domain of lattice sites. For such an operator the general expression of the corresponding 
operators ', aJ f involves M av local av pseudofermion creation and annihilation operators for the av = cO, si 
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branches (and creation operators for the av ^ cO, si branches) which refer to a compact domain of J\f au effective av 
lattice sites. The operators l|46[) - (|53|l involve the product of operators whose expressions involve elementary operators 
of a single av branch and are particular examples of such a general expression, which has the following form, 



and 



Z.LNFJ 



e'v 



-AN. 



II 

3" =3' 



Y\ n fxjii , av fxj,i, , av 

j"=j' j'"=AN" l F + N^ NF +j' 



fx } n , av fxjm . av 

j'"=\ANg F \+N& NW +j' 



av = cO, si 



£d,NF,i 
j',av 



N a „+j'-l 



} I <f x j i / , av ' cv^cO, n<l; sv^sl,m<n. 



Here the spatial coordinates xq, x\,...,XN av -i correspond to the compact domain of Af a v effective av lattice sites 
where the number N av is given in Eq. l(57|) . 

Let us consider the operator Q l ^ au , which is the Fourier transform of the above local operator av , 



This operator can be expressed in terms of the following momentum convolution, 

U k-k',av U av <>k',lAP° v ~ U fc _; Ap0 av U av , 



el- 



where the operator 0jj£'* was given above and 0^^'* is the Fourier transform of the local operator Qj"£f. It reads 



l,NFA 



e 



NF, AT phN F_ 



l,NF,i 
k. av 



= iw) " W a <") n 



i 

Na 



AN"f +N 



i'=0 



n ( e 

<=AN NF +N phNF Qi'=-Qai 



' a av<li' ft 

Jq. 



E « 

AN NF ) 

E 



J q i // , av 



q '' ^i" = AN«F + Ngi NFqt ' 



\AN NF \+Nf hNF -l TV -l 

n n 



+'>:; 



E 



fq^ , av 



■ill - J. 
— %% a av q i n rj 



J Qiii , av 



x S 



> ; av = cO, si 
■>lJ 



p.l,NF,i 
^k,av 



l "^''=° ? * ^." = |AiV < «-P'| + ivgS JVF ' 

i) [n e <-*<**,.- 

x 5j. ; j c ^jv q „-i - j ; cf ^ cO , n < 1 ; ^ si , to > . 



(58) 



Here the set go, <ji, Q2,---,q// au -i refers to j\f av summation canonical-momentum variables associated with the av 
pseudofermion bands. We note that the canonical-momentum values in the Kronecker S's of Eq. (|58(l run under 
the summations but not under the products appearing in that equation. For the av ^ cO, si branches, the J\f a „-av 



pseudofermion operator Q^ N a J l creates and annihilates \AN^ U F \ av pseudofermions when AN^ F > and AN^ F < 
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0, respectively, and generates N^, NF — 0,1,... finite-momentum and finite-energy av pseudofermion particle-hole 
processes. For the av ^ cO, si branches, it creates N^f av pseudofermions of bare-momentum q such that \q\ < q® av . 

In the case of the av ^ cO, si branches, the general expression given in Eq. I|58l) for the operator fc ' is valid 
for densities such that the corresponding ratios n* av — N^ v /N a have finite values. For the cv ^ cO (and sv ^ si) 
pseudofermion branches and electronic density n = 1 (and spin density m = 0) all pseudofermions separate into 2v 
independent —1/2 holons (and 2v independent —1/2 spinons). Therefore, the operator given in Eq. (|58|l does not 
exist. Moreover, for n = 1 (and m = 0) the above generator F^_ h cQ (and F^_ h sl ) of the elementary processes (C) 

reduces to F^_ h c0 — 1 (and F^_ h sl = 1). However, our theory also applies for electronic density n = 1 (and spin 
density m = 0), provided that we consider the corresponding restrictions in the cO (and sl) excitation spectrum and 
take into account that Qif av = ©L^' 1 ) f° r the cv ^ cO (and sv ^ sl) branches. 

For the av = cO, sl branches the pseudofermion weight distribution involves the following matrix element, 

(G\Fj-GS,av Fj-NF,a V F p -h,av ^k] av F G S , av\®) 5 av = cO, sl . 

Our study refers to very large values of L when the commutator [F^_ h au , Fj_ NF av ] = vanishes and the operator 
<=>W is such that 0«£.< F% s = Fl 

GS av an d thus the matrix element can be rewritten as, 

(0\Fj-GS, av F p ~h,au Fj-NF, av ®k-lAP° u ,av F ~GS, a JO) i av = cO, sl . 

When applying the generators Fj l av and F GS av of Eqs. and (JSSJ), respectively, onto the pseudofermion 

vacuum to construct a given energy eigenstate, the set of av band discrete canonical-momentum values {q~j} of the 
pseudofermion operators /| . av in the expressions of these generators are those of the CPHS ensemble subspace which 
that state belongs to. This rule applies when one considers the generators of the full energy eigenstates. (Below we 
express each J-CPHS as a direct product of subspaces. To reach the correct final results, such a rule does not apply to 
some of the states which span such direct-product subspaces.) The same occurs with the discrete canonical- momentum 
values of the pseudofermion creation and annihilation operators of the expression of any operator when it acts onto 
a given energy eigenstate. For instance, let \j3) and |/3') be energy eigenstates. Thus, the discrete canonical-momentum 
values of the pseudofermion creation and annihilation operators of the expression of the operators and O in 0T|/3') 
and ©|/3) are those of the CPHS ensemble subspace which the states \ f¥) and \(3) belong to, respectively. An important 
property for our theory is that for L large both choices lead to the same value for the matrix element (/3|0|/3'). 

The operator 6 fc ' a J l of Eq. is for the av — cO, sl branches such that the commutator [F p ^h,av, a J % ] = 
vanishes when it acts onto the CPHS ensemble subspace which the corresponding excited state belongs to. Further- 
more, there occurs a full overlap of that operator with the generator Fj-jstf, av and for the av ^ cO, sl branches there 
occurs a full overlap of the operator 8t given in the second expression of Eq. (|58|l with the generator FpjF.av 
The latter full overlap results from the lack of av ^ cO, sl pseudofermion occupancy of the initial ground state. 

A J-CPHS ensemble subspace can be expressed as the direct product of subspaces, one for each av branch pseud- 
ofermion occupancy. In the particular case of the av — cO, sl branches the low-energy and high-energy physics 
separate provided that L is large and one can define two of such product subspaces for each branch. They are as- 
sociated with the excitation occupancy configurations generated by the finite-energy elementary processes (A) and 
low-energy elementary processes (B,C), respectively. We call p— h, av = cO, sl branch subspace the latter low-energy 
subspace. Thus, the number of product subspaces equals the number of av branches with finite pseudofermion occu- 
pancy in the J-CPHS ensemble subspace plus two. Finally, for some J-CPHS ensemble subspaces the direct product 
also includes the independent —1/2 holon subspace and independent —1/2 spinon subspace. For such subspaces 
the generator Ft, av either creates N^ v av pseudofermions with limiting bare-momentum values q = iq^v 01 reads 
Fp av = 1 when N^ v — and thus Fj au = F^ NF av (or Fj av does not exist if n = 1 and av = cv or m = and 
av = sv). The states which span the av branch direct-product subspaces and p— h, av = cO, sl branch direct-product 
subspaces have the following form, 

\f.L;av) = Fj_ NFav \GS); av = c0,sl; \f.L;av)=F f NFa JGS); av ^ cO, sl , 
\f.L:p-h,av) = Fl_ htav F\_ GSta JO); av = cO, sl . 

Below we express the spectral functions as convolutions of p — h, av and av pseudofermion spectral functions. How- 
ever, in order to reach the same spectral-weight distributions as by use of the above matrix elements, it turns out that: 

(i) When applying the generators Fj_ NF av and F^ NF av onto the ground state \GS) to construct a given av 
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branch direct-product-subspace state, the set of the av band discrete canonical-momentum values {qj} of the 
pseudofermion ope: 
ensemble subspace; 



pseudofermion operators fq. fCtv in the expressions of these generators must be those of the ground-state CPHS 



(ii) When applying the generator F^_ h au Fj_ GS av onto the pseudofermion vacuum |0) to construct a given 
p — h, av — cO, si branch direct-product-subspace state, the set of the av band discrete canonical-momentum values 
{q~j} of the pseudofermion operators /I. av in the expressions of that generator must be those of the corresponding 
excited energy eigenstate. 

The property (i) ensures that the above full matrix-element overlaps are reproduced. Furthermore, properties (i) 
and (ii) also ensure that the contribution from the unconventional orthogonality catastrophe matrix-element overlap 
discussed below is not counted twice. 

Below we introduce the av — cO, si pseudofermion spectral functions and p — h, av = cO, si pseudofermion 
spectral functions which involve the operators ©^^pF av anc ^ ©a^'* ZAP phF ' respectively. The momentum 
convolution of these two operators leads to the correct expression for the operator ©u* „ such that 8/* = 

1 f±t,NF,i Otl,F,i s _ £*l,NF,i oi.F 

;k lKJ k-lAP F -k'av^av u k'lAP£* F ~ U t-!AP»„ai/ U ai/ 



Sfc' ©fc_iAP ,? -k> av^au' 1 S k' iAP phF = ®fc'-ZAP° au ^w'' In tum ' tlie av ^ c0 ' sl pseudofermion spectral func- 



tions considered below correspond to the operators ^k-uCp" Q „, which refer to the av ^ cO, sl branch direct-product 
subspaces. For the latter branches, the pseudofermion spectral function is associated with the J\f au = > av 

pseudofermions created by the processes (A) , whereas the N F V av pseudofermions of limiting canonical momentum 
iqa V contribute to the independent —1/2 holon (av = cv) or —1/2 spinon (av — sv) spectral function and to 
the momentum of the av = cO, sl spectral functions associated with the elementary processes (A), through the 
FP-scattering-center phase factors. 

The operators associated with the pseudofermion spectral functions can be written in the corresponding direct- 
product subspaces as, 

^k-uLp^au = J2(f- L > av \®k'-iAP F u , a v\ GS ) \f- L > olv)(GS\ ; av = cO, sl , 
/ 

®fc-;AP°„,c«, = J2(f- L > av \®k'-iAPg v , av \ GS ) \f- L "> olv)(GS\ ; av ^ cO, sl , 
/ 

6 k, iap& f = ^(f-L\p~h,av\Q l £' l \GS)5 klApS hF\f.L;p-h,av)(GS\; av = cO, sl . 
/ 

Here the / summations run over the states which span such subspaces and the matrix elements are given by, 

x S . NF phnf_ ; av — cO, sl , 

M av -1 

(f^av\Q l f_^ av \GS) = (^) 2 e-^S-V*, 

x ^ilAP^E^- 1 ^]; ai/^cO, sl, (59) 

for the pseudofermion spectral functions and 

(f.L;p~h,av\Q l ^' l \GS)S k lAp ph F = (0\Fj-GS,au F p -h,au F-GS,aJ°) S k,lAP^ F > av = c0, sl , (60) 

for the p — h, av = cO, sl pseudofermion spectral functions. 

The simple form of the matrix elements (|59|) follows from the full overlap of the generators Fj-np jOCU and Fnf.ov 
with the operator Q^ N a ^} 1 for the av — cO, sl and av ^= cO, sl branches, respectively. Such a full overlap also justifies 
that the corresponding av spectral functions whose expression is given below have a non- interacting character. In 
turn, the evaluation of the matrix element (|60() of the spectral function associated with the av = cO, sl pseudofermion 
elementary processes (B,C) is a more involved problem. For the av = cO, sl branches the phase-factor momentum 
ZAP°„ = l[APP^ F + APjfJ involves a term, IAPP^ F , which arises from the elementary processes (C). Interestingly, 
the dynamics associated with the overlap of the av — cO, sl state (f.L; p — h,av\ = (0\Fj-Gs,au Fp—h, av\ with 
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the state Q^^GS) = F& s> av \0) = F^ GSa jO) of the matrix element (JSDJl is not controlled by the operator 

6.' but rather results from the different discrete canonical-momentum values of the pseudofermion creation and 
annihilation operators involved in the generators of each of these states. (For these branches the expression of 
the operator Q^ av does not include that of the generator F^_ h au , as mentioned above.) Each discrete canonical 
momentum value of the pseudofermion operators involved in the generators of the former state includes an extra 
overall canonical-momentum shift Q a ^(q)/L relative to those of the latter state. If a av — cO, si pseudofermion 
or pseudofermion hole is created at the Fermi points by the elementary processes (B) and thereafter moved from 
there by the elementary processes (C) generated by the operator F p -h,av, the dynamics associated with the overlap 
of the excited-state occupancy configurations generated by the latter processes with the ground-state generator is 
controlled by the orthogonality catastrophe that occurs in the matrix element (|60|) due to the overall phase shift 
Q av {q) /2. Such a matrix element involves N® v + AN F V av pseudofermions. The occupancy configuration of the state 
-GS a JO) corresponds to the densely packed momentum distribution N~®(q~j) for N® v + A.N F V av pseudofermions. 
The corresponding discrete canonical momentum values q~j (occupied and unoccupied) are those of the ground state, 
q~j = qj. The occupancy configuration associated with the state \f.L; p — h,av) — F^_ h aiy Fj_ GS a „|0) also refers to 
N% v + AN F V au pseudofermions. However, its discrete canonical momentum values are those of the excited energy 
eigenstate. This feature leads to an exotic overlap for the matrix element H60I \. Such an overlap is behind the unusual 
quantum- liquid spectral properties, as further discussed below and in Ref. 0] ■ We note that in the absence of the 
av = cO, si pseudofermion overall phase shifts Q av (g)/2, Eq. (|llf> . the matrix element i|60|) would vanish, except for 
the lowest-peak weight such that F p -h,av — 1 and thus APP^ F = 0. 

The spectral function expressions are additive in the contributions of each ground-state - excited-energy-eigenstate 
transition. For each transition, the available excited-energy-eigenstate pseudofermion discrete canonical-momentum 
values are in general slightly different and given by the functional (JSJ) ■ The important point is that for each ground-state 

- excited-state transition one knows the precise values of such discrete pseudofermion canonical momenta. Given these 
values, the av pseudofermion creation and annihilation operators of the matrix element corresponding to the specific 
transition act independently for each av excitation branch. This is behind the introduction of the above subspace 
direct product and follows in part from the factor 8 av> a ' v > on the right-hand side of the pseudofermion anticommutation 
relation l|14l) . Thus, since the pseudofermion creation and annihilation operators of each av branch act independently 
for each ground-state - excited-energy-eigenstate transition, they also do it for the whole spectral function, which 
is additive in the contributions of each ground-state - excited-energy-eigenstate transition. Moreover, as a result of 
the additive character of the energy in terms of av pseudofermion, independent —1/2 holon, and independent —1/2 
spinon single energies and of the corresponding expression of each J-CPHS ensemble subspace as the direct product of 
the above considered subspaces, the excited-energy-eigenstate wave-functions of the ground-state normal-ordered ID 
Hubbard model factorize. It follows that the spectral functions B l ' l (k, u>) of Eq. (|37|l . generated by transitions from 
the ground state to a given J-CPHS ensemble subspace, can be expressed as a convolution of pseudofermion spectral 
functions, one for each branch with finite occupancy in such a subspace and for the independent —1/2 holons and/or 
independent —1/2 spinons, if they have finite occupancy in the same subspace, and two functions for the particular 
case of the av = cO, si branches, as discussed above. It follows from the form of the matrix elements given in Eq. 
(IGUII that the contribution of the corresponding p — h,av — cO, si pseudofermion spectral functions to the weight 
overlaps is more involved than that of the remaining pseudofermion, independent —1/2 holon, and independent —1/2 
spinon spectral functions. 

For each J-CPHS ensemble subspace, we introduce a dimension D associated with the elementary processes (A), 

D = ^9(Jf ca ,), (61) 

where the numbers M av are defined in Eq. I|57l) . For a general J-CPHS subspace the function B l,l (k, oj) of Eq. I|37|) 
can be written as a convolution of the p — h, cO, si pseudofermion spectral function, p — h, si pseudofermion spectral 
function, one M av -av pseudofermion spectral function for each of the D branches such that M av > 0, independent 

— 1/2 holon spectral function, and independent —1/2 spinon spectral function. Thus, let us provide the general 
expressions of the spectral functions corresponding to such a general J-CPHS ensemble subspace. 

The p — h,av = cO, si pseudofermion spectral function associated with the elementary processes (B,C) is given by, 

B Q a A k >") = E \(0\Fj-GS,auF p ^ au Fl GSa jO)\ 2 

J-CPHS-ais-(C) 

x 6{uj-lAEP h /)S kJApSi F ; av = cO, si I = ±1 , i = 0,1,2,..., (62) 

where the energy and momentum spectra are given below and the matrix element is that of Eq. (|60|l . The summation 
^2j-CPHS-av-(C) runs over the J-CPHS ensemble subspace av — cO, si pseudofermion occupancy configurations 
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generated by the elementary processes (C). The indices Q c q and Q s \ remind us that the overall phase shifts of Eq. 
(|llfl have a specific value for each ground-state - excited-energy-eigenstate transition. 

In turn, it follows from the form of the matrix elements of Eq. (|59(l that the av pseudofermion spectral function 
B l a^ F ' l (k,uj) associated with the elementary processes (A) has a non-interacting character and reads, 



B: 



LN FA 



(k,to) 



1 



A'",, 



E 5(w-lAE av )5 k) iAP ai ,; l = ±l, » = 0,1,2,. 

J-CPHS-ais-(A) 



(63) 



both for the av = cO, si and av ^ cO, si branches. Here the summation Ylj-CPHS-av-(A) runs over the J-CPHS 
ensemble subspace av pseudofermion occupancy configurations generated by the elementary processes (A). For the 
av 7^ cO, si branches and densities in the domains < n < 1 and < m < n the number of the latter occupancy 

( Wa jVNi? a ") an( i thus can be written as follows, 



configurations is given by D av 



N a -N + 2£~ =i/+1 (i/ - v)N cv , + 2L C) _ 1/2 - Nl 



Nt 



AfNF 



v>0. 



v > 1 , 



(64) 



where the values of N, N-t, and Ni are those of the corresponding excited-state CPHS ensemble subspace. We recall 
that for the cv ^ cO (and sv =^ si) branches and electronic density n = 1 (and spin density m = 0) the spectral 
function B l c ^ F ^(k,uj) (and B^ F ' i (k,oS)) of Eq. (gSJ does not exist. 

Finally, the form of the operator implies that (GS\Fp av@a F ' l \GS) = 1 for the av ^ cO, si branches with 

finite occupancy in the J-CPHS ensemble subspace. This together with the non-interacting character of the —1/2 
Yang holons and —1/2 HL spinons is behind the form of the independent —1/2 holon (a — c) and independent —1/2 
spinon (a = s) spectral function, which reads, 



• B a!-l/2( fc ' W ) = 7T &{u-lE a ) 8 k ,lP a 



(65) 



Here the coefficient C a is given in Eq. H24fl . While all spectral functions provided in Eqs. i|63[l and l|65|) have a 
non-interacting character, the p — h, cO and p — h, si pseudofermion spectral functions of Eq. I|62l) correspond to a 
more complex problem. The latter functions are further studied in Ref. Q for the metallic phase. 

The av pseudofermion energy spectrum AE av on the right-hand side of Eqs. 116211 and (|63[) can be expressed in 
terms of the bare-momentum distribution function deviations. The energy spectra AE av and AE^f appearing in 
the latter equations and in Eq. H62[). respectively, and the the independent —1/2 holon (a — c) and independent —1/2 
spinon (a — s) energy E a of Eq. Ij65(l read, 

AE av = g AiV:/fe) ea „( 93 ), 



27T 

A£w F = -j- v av [m aUt + i + m aVt _i] 



av = cO, si , 

E a = fi a [L a} _i/2 + S a . + vN aA S a = c, s : 



2^o H. 



(66) 



i/=2 

The momentum spectra corresponding to such energy spectra are given by, 



AP av 



APi 



AP r , 



E AN» F { qj )q j+ AP F v] 

.0 

oo oo 

AJfo + E J - + E J - 



AP£ F = ^[m a „, + x 



m aUt _i] ; av = cO, si , 



Qj=-i: 
Ak F 



v=l 



v=2 



AP S \ = 2k 



aj f -2J2jL 



v=2 



< n< 1, < m <n. 



Pr = 



E AN» F ( qj )[(l + v)ir -qj}; AP SV = E ^/(q^qj; av * cQ, si , 

OO 



8f=-9£ u 
T [-^c, -1/2 



(67) 
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In these equations e cl/ (qj) = 2vfi + e° u (qj) for v > 0, e su {qj) = 2v[i H + e Q su {qj) for v > 1, the bands e av (qj) for 



= cO, si and e°„(gj) for ai/ 7^ cO, si are defined by Eqs. (C.15)-(C.21) of Ref. the small energy AE^ F is such 
that m au , ±1 is the number of elementary <w = cO, si pscudofcrmion particle-hole processes (C) considered above, 
v a u = Vav{q% au ), and v av {q) = de av {q) / dq. 

As further discussed in Ref. |3|, for densities < n < 1 and < m < n the elementary processes (C) leading to the 
spectral- function singular features include contributions from small but finite values of m av . ±i/N a as N a — > 00. For 
n = 1 (and m = 0) the latter processes do not exist for the cO (and si) branch and thus AE^ F = and AP^ lF = 
(and AE^ F = and APf^ = 0). 

Let us consider the general situation when the J-CPHS ensemble subspace has finite occupancy for the cO and si 
pseudofcrmion branches, D — 2 > av 7^ cO, si pseudofcrmion branches, independent —1/2 holons, and independent 
— 1/2 spinons. In this case the functions B l ' l (k, w) of Eq. (|3*7j) can be written as, 

B l '\k,u) = J-G'^w) 



n£)(n(£f"'[ e 

a=c, s a j=l v a/ J-CPHS-oiUj-{A) k' u' 

D D 

x B l £ o (k-l^AP aVj -l Pa-k'^-l^AE^.-l J2 E * 

j — 1 a—c, s j—1 a— c, s 



x B 



where 



^ si (V,c/); G C = C C ; G S = G C C S ; t = 0,1,2,...; l=±l, (68) 



G^(fc,^ = ^^t( fc - fe ^- u 4nFEECft-^4-M 



X 



fel Ul J=l fcj_|_l Wj + i 

-D+2 

[ II B a-,-i/2( fc i - fc j+i> w j- w j+i) JT B q\i ( k D+a,u D+3 ), 

j=D+l a 

°c = ( II C «)[EE r^G^, W )]/[^^ |(/.L;G|e^|GS)| 2 

a=c, s J k f j=l 

Here the coefficient Gc, which also appears in the operator expression of Eq. I|57[l . has a uniquely defined value for each 
CPHS ensemble subspace, 0jy\ j is the corresponding operator on the right-hand side of Eq. (|32|1 . the summation J^f 
runs over all energy eigenstates of the CPHS ensemble subspace, and the summation is over all J-CPHS ensemble 
subspaces of that subspace. Moreover, the pseudofermion spectral functions appearing in the G l ' l (k, u) expression 
are given in Eqs. Ij62|) and (|63|) . the independent —1/2 holon and —1/2 spinon spectral functions are provided in Eq. 
(|65|l . a_D+i = c and cxd+2 = s labels the independent —1/2 holons and independent —1/2 spinons, respectively, the 
momenta ki,k2, kn+3 and energies ^1,0^2, ■■■,ud+3 correspond to summation variables, and the index avj, where 
j = 1, D, is such that avi = cO, av 2 = si, and for j = 3, D avj refers to the D — 2 av 7^ cO, si pseudofermion 
branches such that N^f > for the J-CPHS ensemble subspace. To reach the second expression of Eq. H68|l from 
the expression for G l ' l (k, u>)/Gc, we used the non-interacting form of the spectral functions given in Eqs. il63|) and 
(f 651) to perform D + 2 momentum and energy summations. It follows that the general spectral function B l,i (k, lj) 
of Eq. (|37J) can be written as a convolution of the p — h, cO and p — h, si pseudofermion spectral functions alone, as 
given in the second expression of Eq. I)68[l. We recall that for the i = function B l '°(k, u>) the value of the coefficient 
Gc is independent of U/t and for the dominant CPHS ensemble subspaces considered in Sec. IV corresponding to 
that function, it reads Gc = 1 for all values of U/t. 

For J-CPHS subspaces with no finite pseudofermion occupancy for the av 7^ cO, si pseudofcrmion branches and/or 
no independent —1/2 holon and/or independent —1/2 spinon occupancy, the spectral function B l ' l (k, oj) has the 
same general form as in Eq. JBSJ, except for the absence in the expression of G l,l (k, iS) of the spectral functions 
corresponding to the missing branches and/or quantum object types. Note also that the expression for G l,l (k, oj) and 
thus for B l ' l (k, ui) = G l,l (k, ui)/Gc given in the unnumbered equation after Eq. I|68|) is valid for electronic density 
n = 1 (and spin density m = 0), provided that for j = 3,...,£> the index avj refers to the D — 2 sv ^ si (and 
cv 7^ cO) pseudofermion branches such that N^ F > for the J-CPHS ensemble subspace. Moreover, for n = 1 (and 
m = 0) one must use cO (and si) pseudofermion spectral functions specific to the corresponding excitation spectrum. 
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These functions are studied elsewhere. However, the second expression of Eq. i|t)8|) refers to densities in the domains 
< n < 1 and < m < n only. 

The probability amplitudes |(0|-F_gs, av Fp_ h av^J-GS av\®)\ 2 associated with the matrix element ljfiU|) which ap- 
pear in expression (|62J) have the following general form, 

2 



\F-GS, a uF, 



p—h, olv J—GS 



F 



,aJ0> = (0|/ g> 



xv " " " fq' j , av fq y ,av ' ' ' fq 



Q N o +4s f ,a 



JO) 



(69) 



where av — cO, si. In expression (|69|) we have considered that the ground state corresponds to pseudofermion 
annihilation operators and the pseudofermion operators left for the excited energy eigenstate are of creation character. 
At this stage, for the evaluation of ./V-electron spectral functions, the main problem remaining is the computation of 
the non-trivial probability amplitude (j69JI . which can be expressed by the following determinant, 



{fqi,av t fq\,av} 



{f\i,av i fq' x ,av\ 
{4„o . »»p ,otvi fq\,av) {fq 



{f\l,av i fq' 2 .av} 
{/§2> av ) fq' 2* av \ 



{/ 8i, av ■ J'i N v \ v 

{fq2 , OLV 1 f<H' N 



av\ 

F ) OCVj 



v ) fq' 2 , av \ ' ' ' {fq. 



,av > fq' N o / + A«f / ' t "'} 



(70) 



for the av = cO, si branches. This result is justified by the following pseudofermion properties. First, the pseud- 
ofermions have no residual- interaction energy terms, as discussed in Refs. tjj. |!|. Second, the canonical-momentum 
shift Qavilj)/ L, which under the ground-state - excited-energy-eigenstate transition involves all av = cO, si pseud- 
ofermions of the initial ground state, is a zero-energy process |9J. Third, the elementary pseudofermion processes (C) 
correspond to av = cO, si pseudoparticle particle-hole processes whose energy spectrum is of non-interacting char- 
acter for the pseudoparticles, what implies that the energy spectrum AE*^ = [27r/L] v av m av of Eq. (|66|) remains 
linear in m a „ for small finite values of m av /N a as N a — > oo |3j- 

In spite of the non-interacting form of the determinant l|7U|l , the unusual pseudofermion anticommutation relations 
(|14|l give rise to unusual physics, in the form of an orthogonality catastrophe. (The absence of such an orthogonality 
catastrophe would require that Qav(o)/^ = 0.) Indeed, replacement of the anticommutator i|14fl in the determinant 
(t7T)|l leads to, 



sin y °2 — — J 

7 92— ¥ 
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y_^u _jut I Bin^ I sm^^" 



(71) 



for the av = cO, si branches, where the overall phase shift Q a v(qj)/2 is given in Eq. (|llfl and the bare-momentum dis- 
tribution function N"^,{qj) is such that N"^,{qj) — Mf^,(qj). Here Navies) ^ s the pseudofermion canonical-momentum 
distribution function given in Eq. 156(1 . which includes the low-energy and small-momentum av — cO, si pseud- 
ofermion particle- hole processes (C). The determinant of Eq. H71fl can be rewritten as, 
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(72) 
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where N~®(qj) = Af^{q_j) is a densely packed bare- momentum distribution function whose Fermi points are given 
by qFav, ±1 = ±<7fc<j; ^ [2i'/-^]^^'olf±i an d the corresponding canonical-momentum distribution function N"av(<lj) * s 
that of Eq. (jS7j)l whose Fermi points read gpw,±i = i^Fc^ ^ [^/^[AA^,, ±1 ± QtLC^ffFai/V^ 71 "]- The expressions 
CHJ-C2J are used in Ref. Q in the derivation of finite-energy spectral-weight distributions for the model metallic 
phase. 

VI. DISCUSSION AND CONCLUDING REMARKS 

The main result of this paper is the general spectral function expression defined by Eqs. H36|) and (|68[) . The expres- 
sion given in the latter equation involves the p — h, cO and p — h, si pseudofermion spectral functions provided in Eq. 
iJB^J, whose probability amplitude | (0| Fj^qs, au F p -h, au F~gs aJ0)l 2 can ^ e expressed in terms of the determinants 
of Eqs. H70|l - H72|) . An important aspect of the pseudofermion dynamical theory introduced in this paper and further 
developed in Ref. 3] for the metallic phase, is the different origin of the dynamics associated with the matrix-element 
overlaps of the au = cO, si pseudofermion occupancy configurations in the vicinity and away of the Fermi points. 

In reference [2| the finite-energy spectral function expressions derived by use of the pseudofermion dynamical 
theory introduced here are applied to the study of the spectral- weight features observed in the quasi- ID organic 
compound TTF-TCNQ. Interestingly, one finds quantitative agreement with the observed spectral features for the 
whole experimental energy band width. The microscopic mechanisms found in Ref. by use of our theory are also 
consistent with the phase diagram observed in the (TMTTF) 2 X and (TMTSF) 2 X series of organic compounds and 
explain the absence of superconducting phases in TTF-TCNQ. Our theory is also of interest for the understanding^ 
the spectral properties of the new quantum systems described by ultra-cold fermionic atoms on an optical lattice |(| . 
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